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i UCTIO 


By the term “Remote Sensing" we mean the procedure by 
which one can get information about a target and its 
environment with the use of a sensor which is not actually in 
contact with the object under study. 

In our modern era one of the most important methods for 
surveillance and tracking of targets is infrared (IR) 
technology. Infrared technology consists of devices whose 
information depends upon the electromagnetic radiation of 
wavelengths between 0.7 and 1000 micrometers (0.7 to 1000x10- 
4 cm) which is absorbed by, reflected from, or emitted by 
different objects. By using the natural emission from objects 
most infrared systems can be used as detectors or target 
designators. Since infrared systems are passive the 
advantages of such systems are enormous. 

But as for all systems the infrared systems have their 
limitations. Particle effects on transmission through clouds 
and fog limiting the performance of the infrared radiation 
almost as much as the visible radiation, although the 
infrared can penetrate the atmospheric haze more effectively 
than the visible mene because of its longer wavelength. The 
various atmospheric gases absorb at several wavelengths of 
the infrared region leaving some areas where radiation can be 


effectively transmitted, known as atmospheric windows. 
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Since all targets which have a temperature above absolute 
zero radiate energy there is always a certain amount of 
radiation from the background present that must be contended 
with for the detection of a weak or a distant target. The 
reflected sunlight in the near infrared region (< 3 
micrometers) represents an important problem in the detection 
of a target (aircraft) against a cloud background at daytime. 

The above factors imply that a detection and surveillance 
infrared system will have a limited all-weather or daytime 
capability unless some techniques are brought to bear which 
will discriminate the target from the background. 

For the development of such a system it is important to 
specify the infrared radiance for an earth’s atmospheric 
scene. What is also required is the full knowledge of the 
spatial and spectral structure of the background as well as 
the atmospheric effects on the emitted thermal radiation and 
intensity of the target. 

The radiation reflected or emitted from the target mixes 
with the radiation from the surrounding objects and creates 
confusion for a detection and surveillance system. 
surrounding objects may include clouds, buildings, land etc 
which reflect radiation from the sun or emit radiation of 
their own. The sea radiance which is largely due to the 
sun/sky radiation reflected by the sea must also be taken 
into consideration in designing a system for target detection 


for low angles above the sea water (ships or low flying 
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airplanes). Most systems for detection and recognition depend 
on the "Radiance Contrast" described by the "Temperature 
Difference’ between the target and the background. It must be 
noted that the "Temperature Difference” can be deceptive, in 
the case of the radiance contrast, unless the reflected 
radiation is also considered. 

For thermal imaging systems the background scene’s 
radiance is usually suppressed. That suppression allows the 
creation of high image contrasts. 

A good detection and surveillance system can be created 
only if magnitudes and variations in the radiance from 
different background scenes have been studied and taken into 
account with the infrared target appearing in them. 

At the present time NACIT is developing an Infrared 
search and Target Designation system (IRSTD) modified from 
the original ADM version of the AN/SAR-8 (TRS®) which®@will 
provide threat warning and target track information. One of 
the potential problems for the IRSTD is the false alarms that 
may result from the presence of clouds and sea clutter in the 
environment of the target. It is important to understand that 
an infrared sensing system will try to detect the infrared 
radiance in a specific field of view (FOV) requiring high 
spatial resolution data. For this reason NSWC and Naval 
Research Laboratories (NRL) in support of the Naval 
requirements for fleet defense systems created a computer 


program for an IBM/AT micro computer called the NSWC Clutter 
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Measurements Code which was further supported by the 
Background Measurement and Analysis Program (BMAP) for the 
analysis of background scenes appropriate to generic infrared 
systems. A copy of that code, generated by ONTAR Corporation 
under contract, was made available to NACIT. This code 
includes three stored measurement data sets, other processing 
techniques, a version of the LOWTRAN propagation code and a 
high resolution graphics package. The program uses new 
improved techniques in signal processing so that background 
effects may be minimized and an IR surveillance/detection 
system may improve its eaneb i eae: and optimize its 
performance. 

The purpose of this thesis is to perform an evaluation of 
the BMAP code to provide future data to the IRSTD data 
processing and sequence aaton system which is now in its final 
Stage of engineering development. For that reason the 
adaptation/use of the BMAP program was made by transferring 
data from the main program to a Masscomp computer associated 
with the IRSTD. The image data were located in three files 
called BEDFORD files, named after the place in which the 
experiment took place, Bedford Ma. The data was collected 
there from 10 to 14 September 1984. The observations were 
made at a hole between clouds in the sky backlighted by the 
Sun using a dual scanning radiometer in the 4 to 5 micrometer 
and 8 to 11 micrometer bands. The detectors consisted of 16 


element arrays. For the experiment it was assumed that the 
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clouds were constituted of water droplets with a logarithmic 
Size distribution and a mean radius of 4 micrometers. 

The present work consists of ten chapters and four 
appendices after the present introduction. In the first and 
second chapter a theoretical overview of the background 
problem is presented. In the third chapter the problems, 
related to target .tracking, of the various background 
radiances are discussed including references to the marine 
background. The fourth and fifth chapter are about the main 
NSWC clutter code (in conjunction with the BMAP source code), 
presenting an adaptation of the program with its full 
capabilities including the LOWTRAN6 code, and highlighting 
some important aspects for the future user. In the sixth 
chapter a discussion of the collected data with the BMAP is 
presented including the method to enter that data into the 
Masscomp part of the IRSTD using different computer 
algorithms. The sixth chapter presents the format in which 
the image data are arranged inside the Bedford files so that 
they can easily be transferred from one system to the other; 
also a method for extraction of the image data from the 
BEDFORD files is discussed. In the seventh chapter we have a 
presentation of the organization of data inside the Masscomp 
computer of the IRSTD including the method of transfer of the 
data, already extracted from the BEDFORD files, to magnetic 
tape in the Masscomp computer. In the eighth chapter a 


connection between the data format of the portable AGA 
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Infrared System and the formats of the NSWC clutter code and 
MASSCOMP is made. In the ninth chapter an analysis of the 
method of extracting the data from the BEDFORD files is made. 
Finally in the last chapter conclusions and recommendations 
for future work are made. 

In Appendix A the LOWIN and LOWPLOT cards of the 
mainframe program LOWTRAN6 found as an option inside the BMAP 
program are presented; in Appendix B the computer algorithm 
to extract the image data from the BEDFORD files is listed. 

Appendix C provides a useful description ina few words, 
for the future user, of the “"F1...F10" keys of the IBM 
computer in use with the NSWC clutter code program. Finally 
Appendix D contains the computer software initially provided 
by the ONTAR company to read the BEDFORD files. 

It is expected that the results from this work will 
provide a start to further development of the BMAP program 
for use with the IRSTD system (always considering the limited 
capabilities of a small portable computer) by transferring 
data from other infrared detection systems, such as the AGA 
Thermovision 780, into the BMAP and later using the computer 
routines of this thesis to transfer the infrared image data 
to the IRSTD. Furthermore it is expected that this work will 
@ive a good understanding of the background problems which 


may appear in an environmental infrared scene. 
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A. BACKGROUNDS AND TARGETS 

An infrared system must receive the radiation from the 
target, identify it as aninput signal, and distinguish it 
against all background radiation within the field of view. 


But since all objects which have a temperature above absolute 


zero radiate energy in the infrared, that background 
radiation iS omnipresent, causing noise in the optical 
system. This is called the external noise and has two 


different components [{Ref.3]. The first component is called 
the spatial component, caused by the different structural 
details of the environment and the second component is called 
the temporal component, which is actually a photon noise, 
caused by the fluctuations of the photons emitted from the 
background (proportional to the square root of the number of 
the photons). These photons cause a modulated flux on the 
detector called "“Background-Radiation Photon Noise”. The 
lower limit for the detector’s noise, and hence the best 
detectivity, is determined by these fluctuations. Using the 
Planck’s law and the Bose-Einstein statistics, knowing the 
environmental conditions, we can calculate the background 


photon noise. 


B. SPATIAL BACKGROUND RADIANCE (NOISE) 
The spatial background noise can be produced by many 
factors such as sunlight scattered from atmospheric particles 


ed 


like water droplets, fog, haze and clouds, or from sunlight 
reflected from the sea surface/terrain or self emission of 
different objects like the sun, stars and planets or objects 
on the surface of the earth like people, trees, etc. There 
are two different kinds of spectral background radiance: 

i Celestial background 

ae Terrestrial background 

The celestial background is the radiance due to the sun, 

moon, and planets which affect space-based optical systems 
looking out to space. This is important for many applications 
but not for the subject of this thesis. 


The terrestrial background is the radiance coming from 


various sources inside the earth’s atmosphere. This 
background can be considered to arise mainly from three 
sources: 

le the scattering of solar radiation by the earth or 


clouds or the sea, 

Ze the emission from different molecules in the 
atmosphere such as carbon dioxide, ozone molecules 
and water vapor, and 

oe the emission of surface objects like trees, power 
plants, sand, people, houses, etc. 

The full analysis of various terrestrial backgrounds is given 


on Chapter III of this thesis. 
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C. TARGET-BACKGROUND BRIGHTNESS/TEMPERATURE DIFFERENCE 

One Cf Mene omost important factors in the detection 
probability of a target is the radiance difference between a 
target and its surrounding background. That radiance can be 
expressed with the radiation contrast described in Ref.3 with 


respect to the brightness as: 








Bo - Bb 
Can tm when Bo > Bb Go-4.) 
Bp 
or Bp - Ee 58 
CO BS_.. = ee when Bo < Bo (2-2) 
Bo Bb 
where C = radiation contrast 
Bn = average brightness of the background where the 


target is located, 

Bo = target’s brightness, and 

5B = incremental brightness of the target from the 
average background Bb. 

The contrast C may vary from 0% to 100% if the target is 
darker than the background. If, on the other hand, the target 
is brighter than the background then the C may vary from U to 
Paranaty. 

The detection of a target against the background depends 
on the background brightness and the size of the target. 
Figure 2.1 [{Ref.3] presents the threshold contrast vs. 
background brightness for various angular sizes (in minutes) 


of the target. 


oH 


Another method of describing the contrast factor [Ref.1], 


named the "Radiation Contrast Factor” Cr, is in terms of the - ~— 


apparent temperature difference between the target and its 


background environment. Cr is defined from the equation: 


Wr - Ws 
CR -=-——— (Zo) 
Wr + WB 
where Wap = Background radiant emittance (W/m2) and 


Wr = Target emittance (W/m2). 
Since Wr=¢roTr‘4 then: 
CR = 265T/TB foreol<<T (2-4) 
where 65T = Wr - Ws 


Tr = Target temperature 


Ta = Background temperature 
3 
Me Target Size (minutes) 
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Figure 2.1 Threshold Contrast versus Background Brightness 
(From Ref.3) 
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In the case of small temperature differences between the 
background and the target, the radiation contrast factor 
becomes very small and the detection of the target becomes 
weny difficult. 

The radiation contrast as a function of the target-to- 
background temperature difference is shown in Figures 2.2 and 
2.3 describing the radiation contrast for the 3 to 5 and the 
8 to 14 micrometers infrared band. In these figures the 
vertical axis represents the Radiation Contrast factor and 
the horizontal axis represents the temperature difference 
between the target and the background for various background 
temperature values (from 240 K to 300 K). 

A good example of radiation contrast is presented by 
Seyrafi in Ref.3 describing an experiment involving a truck, 
parked in an open field and observed over a 24 hour period 
With an infrared optical system. In the experiment it was 
observed that the contrast between the vehicle and the 
background varied from positive values in the daytime to 
negative values during the night passing through zero (equal 
temperature) twice a day. More specifically the truck in the 
afternoon gave a positive contrast value which increased up 
to a maximum — At night, after passing from the equal 
temperature value, the truck was cooled faster than 
background so that the contrast value became positive again 


and the cycle continued. As the difference in temperature 


23 


between the target and the background increases the potential 
between the target and the background increases the potential 
acquisition range increases too as presented in Figure 2.4 
[Ref.4]. In the figure we show the acquisition range for a 
target as a function of target-to-background temperature 
difference. Also the increase in the temperature difference 

leads to an increase in the range for identification and 
Classification of a potential target. In Figure 2.5 [{Ref.4] 
we have a range prediction for identification and 
classification of a target with respect to target-background 


temperature difference in degrees Kelvin. 


D. SPECTRAL-SPATIAL DISCRIMINATION (CLUTTER SUPPRESSION) 

The process of encoding the target’s information against 
the background suffers from the problem that the signals from 
the background will also be processed and modulated like a 
real target. A system is said to be background noise limited 
when the unwanted signals or noise from the environment 
exceed the other noise sources of the system. When we scan 
the different background objects, the noise created will 
depend on their spatial distribution, spectral radiance, and 
on the encoding and processing method employed. 

In one target discrimination model [Ref.5] a distribution 
of background noise amplitude is obtained and a threshold 
noise level set on the basis of expected frequency of 


occurrence according to environmental conditions set. 
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Figure 2.2 Radiation contrast for the 3.5 to 5 um band 


(From Ref.3) 
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Figure 2.3 Radiation contrast for the 8 to 14 um band 


(From Ref.3) 
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Figure 2.4 Target Acquisition Range as a Function of Target 


to Background Temperature Difference (From Ref.4) 
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Figure 2.5 Target Classification and Indentification Ranges 


as a Function of Targey to Background Temperature 


Difference (From Ref. 4) 
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Another model [Ref.5] uses the “Wiener Spectrum" which is a 
Fourier transform of a two-dimensional autocorrelation 
function. This is the most effective method of describing the 
spatial variations of background radiance with respect to 
confirming the mathematical model with the real physical 
conditions. 

Finally we can say that the value of the limita noise 
is largely dependent on the electronic filtering, the spatial 
filtering, the spectral filtering and all the other filtering 
that. «can be applied to the signal regardless of the 
mathematical model used to describe the background noise. 
That is the only way to achieve a background noise limited 
system. 

In infrared devices we have a spectral discrimination 
which is due to the spectral response of the detector/filter. 
To achieve maximum average contrast we use different 
statistical methods depending on the application. If we 
attempt to detect a target signal S(%3%)={t(3%) - b())] where 
t(X) is the target emitted power in the presence of a 
background random signal b(12A) and taking in account the 
preamplifiers noise N (independent of the choice of the 
filter) then to increase the average contrast between the 
target signal and the average background signal we must have 
a filter whose function f(2%) must maximize the ratio in 


equation (2-5): 
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A filed eeoly ) dd 
Cao) 


> Ope OX) dA) + N 


Where S( A) is the average radiant power as a function of 
the wavelength, b(.\) is the average variation in background 
radiant power as a function of wavelength. (A complete 
presentation of this filter analysis is presented by Khalil 


Seyrafi in Ref.3). 


E. WIENER SPECTRUM/POWER SPECTRAL DENSITY 

The power spectral density function, often called the 
Wiener spectrum, is a representation of the statistical 
characteristics in the frequency domain of the random input 
autocorrelation function. In a time-invariant system, for 
linear response ina stationary random process, we usually 
use the power spectral density function to derive the mean 
square response for the input variable and the appropriate 
transfer function. 

In the specific case of a tracking loop analysis the most 
practical method for-time response is to represent the inputs 
by random time series input signals and to use Monte Carlo 
statistical methods, as will be presented later, to find the 
power spectral density function. 

We will specify the average domain with "w" [Ref.5] and 


we identify a sample Sw(y) in a random scene radiation 


eo 


distribution Sw(x) for a two dimensional process. The average 
@(x) (an autocorrelation with y a displacement variable in 


the same direction as x) can be written as: 


+ © 


G(x) =-@ ow(y) Sw(xty) dy (2563) 


The fourier transform of $(x) is called the Wiener 


Spectrum of the scene and is defined with W(k) as: 


+o 
Wik) ==! G(x) e-mik.x dx (2-7) 


The Wiener spectrum actually gives a full description of the 
Statistical properties of.a gaussian process and it isa 
function of two variables but can be specified completely as 
a function of a single spatial frequency variable. 

The autocorrelation function and the Wiener spectrum are 
a fourier transform pair because both contain the same basic 
information about a random process, but in practice the 
Wiener spectrum is a easier representation since the 
convolution of the space functions is replaced by a product 
of functions of spatial frequency. So if we want to analyze a 
system’s performance:it is better to use the Wiener spectrum 
in order to manipulate aperture and image functions than to 
use convolution equations. 

The Wiener spectrum can be measured in many ways. Figure 


2.6 {Ref.7] provides a representation of various methods of 
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Wiener spectrum 
measurement 
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Figure 2.6 Techniques for measuring Wiener spectra 


(From Ref.7) 


measuring this spectrum which can be implemented either by 
digital or analog techniques. 

To get a significant target signal its surrounding 
background we must maximize the ratio of the instantaneous 
target signal squared divided by the mean squared background 
Signal which is the following expression [Ref. 8]: 


+c 2 
| va A*(k) T(k) dk 





(2-8) 
ee 
- © | A(k)|2 We(k) dk 
where T(k) = Fourier transform of target’s distribution, 
A(k) = Fourier transform of the detecting 


aperture, and 


Wap (k)= Wiener spectrum of the background. 
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By using the background Wiener spectrum and the target 
radiance function we can create various space filtering cases 
in which the spectrum must be selected according to the 
environment Ou the target and Similar statistiea: 
characteristics. We are actually concerned with the peak- 
Signal to rms background-noise ratio at the output of the 
system radiation detector. For that reason we must choose a 
mathematical description of a space filtering field that will 
give us the wanted maximum. 

The Wiener spectrum has some limitations too. The first 
one is the assumption that the spectrum has a quadratic 
probability density. Under that assumption we come up against 
the problem that a quadratic function for some negative 
values of the argument must take real values. But we can’t 
give negative values to the spectrum of the vaaiance because 
it is always positive valued. The best we can do in solving 
the problem is to consider zero values in the place of the 
negative values of the argument for the random radiance 
process. 

The second limitation is that when we examine the 
instantaneous field of view (IFOV) of a system for a cloudy 
environment we cannot apply the central limit theorem [Ref.8] 
because the detector can capture only a few shots of the 
cloud form and the characteristics of the infrared background 
fluctuations. As a result we cannot determine the higher 


moments of the joint probability functions. We can however 
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make better predictions of peaks of the background noise by 
measurements of the statistics of phase relationships and 
higher moments. 

We can take for example a partially clouded sky and prove 
that a background scene must be described by more than one 
principal peak in the power density spectra. The scene 
considered could for example have a bimodal distribution with 
two peaks in which one at the lower radiance values describes 
the space between the clouds and the other describes the 
scattering and reflection in clouds with a range of radiance 
level of high probability. In that situation we could have 
problems describing the probability density function with a 
gaussian curve because it would not represent the actual 
scene.’ These problems are most obvious in the case of 
calculations of false probabilities in which the enee 


spectra are not appropriate background descriptions. 


F. MONTE CARLO SCATTERING 

The Monte Carlo procedure is a statistical method to 
describe a complicated set of physical processes of photon 
transmittance in the atmosphere. 

When the Bear and the emitting source directions are 
determined we need a large number of trial photons because 
until the calculation is finished we do not have any 
information about the exit direction of each photon. We can 


improve this by accumulating the exiting photons in angular 
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bins. We can track each photon until it makes a predetermined 
number of scattering events or leaves the cloud. The ratio of 
the scattering to the total cross section (IR single scatter 
albedo) is smaller than 1.0 (the ratio usually varies from 
0.5 to 0.9) so that after a finite number of scattering 
events (about 20 to 50), the photon becomes insignificant and 
can be dropped. 

To describe Monte Carlo scattering we <u select random 
values in the range (a,b) which are then normalized in order 
to be inside the panes O to 1. That can be done by using the 


equation: 


+o 
a i(a,*%) ee ae 
a C299 


+o 
a) 2 (a; Oe (so aes 


Function defining the values of the parameter, 


where F(x) 
and B = randomly selected value. 

We then invert the equation to give the corresponding value 
OL 2G. 

The whole process must be made in a reasonable amount of 
time to obtain meaningful results and is described in the 
following way: 

gi We first calculate the spectral weighing function 
from the atmospheric transmittance, solar irradiance 
and instrument response. By random selection in the 


interval we calculate the wavelength of each photon 
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and then from the tabulated values we determine its 
albedo and its scattering phase function. 

The photon propagating through the cloud for a 
distance R, determined by a Monte Carlo selection, 
encounters an aerosol. We then make a check to 
investigate whether the potential scattering event is 
inside the cloud. We continue the calculations from 
step 4 if we find that the photon has escaped from 
EMcwc1OuUGad, cither from the Sides, top, or bottom. . 
From the two Monte Carlo selections Pore. tne 
scattering angle we determine a new direction for the 
scattered photon. In the first Monte Carlo selection 
we pick the azimuthal angle with the assumption that 
all values are equally probable. In the second Monte 
Carlo selection we pick the scattering angle using 
the Henyey-Greenstein function (Ref 5). We then 
repeat step 2 with the new direction from the 
Seaucering DOINnt . 

When the photon exits the cloud it is binned 
according to its escape angle. We then tabulate the 
reflected and transmitted escapes and obtain the 
total scattering probability by assigning a weighting 
function Wn in which W is the single scattering 
albedo and the n is the total number of scattering. 
The procedure is then repeated for a new photon and, 


when we have completed the process for the assigned 


aD 


number of photons, we calculate the distribution 
functions. 

6. We can predict the amount of the incident energy 
coming out ina certain direction by the number of 
the photons in the same angular bin. We then select 
the bin and get the scattered intensity for the 
specified observer direction and as a by product we 
also get the distributions for all the other 
directions too. We can also improve the statistical 
uncertainty, when it is too big, by averaging the 
distributions over several angular bins. 

The method by which the Monte Carlo method is implemented in 


the BMAP program is presented in Chapter V of this thesis. 
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Til. RADIANCE IN THE BACKGROUND PROBLEM 


A. SPECTRAL RADIANCE OF THE SKY 
1. Introduction 
For a good analytical presentation and understanding 
of the background radiation problem it is important to have a 
good knowledge of the sources that cause it. Sky radiance is 
caused mainly by two sources: 
a) scattering of solar radiation, and 
b) emission of the various atmospheric molecules such as 
CO2, ozone and water vapor. 
2. Scattering 
The scattering of radiation in the atmosphere occurs 
by extracting energy from the incident light and, taking the 
center of the particle as origin, reradiating that energy 


into a 4x solid angle. 
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Figure 3.1 Scattering patterns (From Ref.6) 
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This scattering depends on the particle’s isotropy, 


refractive index and its size relative to the wavelength. In 


Figure 


3.1 {Ref.6] we have three scattering patterns from 


three different sizes of particle. 


a) 


b) 


Co} 


When the particle is smaller than 1/10 of the 
radiation wavelength, we have symmetric 
scattering in the backward and forward direction. 
When the particle is of size 1/4 of the lights 
wavelength we have larger scattering in the 
forward direction. 

When the particle is of larger size than the 


wavelength we start observing side lobes. 


The wave which is scattered is polarized by a small amount 


but since the wavelength remains unchanged we have no loss of 


energy. The scattering phenomena can vary according to the 
size of the particles found in the atmosphere which are: 
ones RADIUS (um) CONCENTRATION (No/cm? } 

Air molecule 10-4 1019 

Aitkem nucleus O22 10-7 104-102 

Haze particles £_£$=10-2-1 G0 

Fog droplet ae 10G@=10 

Cloud droplet r=9 30 0R 0 

Raindrops 102 =103 1072.0 2 
The case of molecular scattering is most commonly known as 


the “Rayleigh scattering" and varies proportionally to the 
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second power of the particle volume and is inversely 
proportional to the fourth power of the wavelength [Ref.10]}. 

The other case of scattering is that caused by larger 
particles (usually of irregular shape) and is usually known 
as “Mie scattering"; here the common assumption is that all 
atmospheric particles regardless of size may be considered as 
isotropic spheres. In the limit of small sizes Mie scattering 
includes Rayleigh scattering. 

o. ion o t £3 

According to Kirchoff’s law [Ref.6] € = I/Iv where € 

is the radiant emissivity of a particular surface, [I is the 
total radiant power per unit area emitted from the surface 
and Ib is the total radiant power per unit area emitted from 
a blackbody of the same temperature. The atmospheric gases of 
various species produce emission of radiation according to 
their absorption bands. For the most common temperature found 
in the atmosphere almost all of the radiant energy is between 
3 and 30 um with the maximum (blackbody) peak occurring at 10 
um With a value of 107-3 watt cm-2 um! steradian!. The 
reason is that the effective temperature is usually at the 
range from 200 to 300 K, so that the blackbody radiance gives 
the most probable value for the emitted radiance. 

The maximum radiance of the sun occurs at a wavelength 
of about 0.5 um, since it radiates like a 6000 K blackbody, 


With a value of 3x1i0-2 watt u-! steradian-! [{Ref.6]. 
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Figure 3.2 Spectral radiance of the sky in a idealized form 


(From Ref.9) 
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[It has been found experimentally that the scattered sunlight 
radiance and the molecular emission radiance in the range 3 
to 4 um are equal [Ref.6]. 

When we go to shorter wavelengths the molecular emission 
can be ignored in the daytime and has avery small value at 
night because the scattered sunlight has a much larger value. 
The situation changes when we go to wavelengths above 4 um. 
There the scattered sunlight can be ignored since the 
molecular emission is much greater. In Figure 3.2 [Ref.9] we 
have a presentation of the spectral radiance of the sky with 
the sun radiating,as mentioned before, as a 6000 K blackbody 
and the sky as a 300 K blackbody. From the figure we can see 
that we have a uniform diffusion of the radiance due to solar 
radiation presented by the "sunlit cloud curve”. 

The average radiance is approximately 2x10-5 times the 
sun’s radiance, although practically only 1/10 of the maximum 
spectral radiance is scattered. In the figure this is 


presented as the "clear sky curve’. 


B. SPECTRAL RADIANCE ABOVE 3 um - CLEAR SKY 

The spectral radiance is a product of a blackbody’s 
radiance, varying with temperature, and of the emissivity 
depending on the number of the emitting molecules along a 
path. Hence the spectral radiance of the sky strongly depends 
on the following: 


1) air temperature 
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2 ) water-vapor content of the atmosphere 
is} ozone, and-~--- 
4) Steyvanion angle from the horizon. 

The dependance of sky spectral radiance on elevation 
angle is presented in Figure 3.3 [Ref.6]. In that figure we 
have the spectral radiance of a clear night sky at 11000 feet 
above sea level and ambient temperature of 8° C for 
elevations of: 0°, 1.89, 3.6°, 7.2°, 14.50, 30° and 90°. From 


the figure we can see that we have only weak absorption lines 


in the 8 to 13 um region so we have some slight dips. 
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Figure 3.3 Spectral radiance of a clear night sky at Colorado 
(From Ref.6) 

For the path at 0° elevation we have a finite air mass 

because of the earth’s curvature and we have almost blackbody 


radiation (€=1). If we increase the elevation angle we see 
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that the spectral radiance decreases and we no longer have 
blackbody emissivity. The emissivity for elevation angles 
near the zenith is low at the regions from 8 to 13 um, and 
the spectral emissivity is almost unity for the water vapor 
band of 6.3 um and the carbon dioxide band of 15 um. We can 
also see a peak emission at 9.6 um because of the ozone. 

The spectral radiance is highly dependent on the ambient 
temperature. From Figure 3.4 [Ref.6] we can see the spectral 
radiance of a clear sky with an ambient temperature of 27° C 
and elevation angles of 0°, 1.8°, 3.6°, 7.29, 14.50, 30°, and 


90°. 
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Figure 3.4 Spectral radiance of a clear sky at Florida 


(From Ref.6) 
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Comparing this 


figure with Figure 3.3 we can tell that the 


humid atmosphere of Florida at sea level makes the emissivity 


go to 1 for the wavelengths less than 7 um and higher than 15 


WM, 


of water vapor and C0O2, 


because we have a greater air mass with larger quantities 


but in the middle region of 8-13 um 


we have a higher emissivity for all elevation angles. 


In Figure 3.5 we have a comparison of the spectral 
radiance for different ambient temperatures at Colorado 
(2. 5om€ and 27 25emG): 

The dashed lines represent the black body spectral 
radiance for the ambient temperatures at the time of 
measurement. 
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We can see that we have spectral radiance almost at the 
blackbody limit for the wavelengths corresponding to the 
absorbing bands and completely different when we are out of 
this region because of the smaller emissivity and colder air 
mass since the atmospheric layer is at a greater distance 


from the observation point. 


C. CLOUD RADIANCE 
The most common objects in the earth’s atmosphere are 
clouds. The clouds can be separated according to altitude 


into three major types with subsequent variation as follows: 





LEVEL HEIGHT OF EXISTENCE Aver OMmeLOUDS 

Low Earth’s atmosphere to 2 Km Stratus-Stratocumulus 

Middle ae i) «(Km Altocumulus 
Cirrostratus 

High em 13 Km Cirrocumulus 
Cirrus 


Their behavior is governed by the following rules: 

1) they can scatter solar radiation, 

2.) they can emit radiation from their own body due to 
larger temperature from their surrounding 
environment, and 

5 they can reflect radiation which originates from 


terrestrial objects. 
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Figure 3.6 Thermal Radiation through clouds (From Ref.11) 
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ieee Clouds contribute a major portion of the spectral 
radiation of the sky since by scattering the sunlight they 
become a source of radiation themselves. In order to do this 

they have some limitations of their own. They must not beless 
than 100 um in thickness and they must have droplet densities 
less than 0.5 g/cm2; otherwise they behave like a blackbody 
with emissivity and absorption coefficient approaching unity. 

Figure 3.6 [Ref.11] shows a schematic presentation of the 
reflection of thermal radiation by clouds. About 20% is 
reflected and 3% is absorbed while 26% of the radiation going 
back to space is originated by the clouds. 

In Figure 3.7 [Ref.6] we have the spectral radiance of a 
dark cumulus cloud in which the dashed lines represent the 
cloud temperature (-10°C) and the ambient temperature of a 
blackbody. 

When we come to higher elevations (cirrus clouds) the 
spectral radiance is represented by Figure 3.8 [Ref.6] for 


various elevation angles. 


D. RADIATION OF AN OVERCAST SKY 

The mieuds: in .the atmosphere are composed of water 
droplets which absorb, refract and reflect the infrared 
radiation. In the visible region the clouds form a good 
diffuse reflector because the water droplets do not absorb 


any radiation. Figure 3.9 [{Ref.1i]}] shows the spectral 
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emissivity vs. wavelength for cloud layers of thickness 6z 


— =-- 
‘oe 


~- 
—e eee - 


measured “in meters. 


When we have an overcast sky, in which the highest cloud 


is at an elevation of 500 to 700 feet, then the radiance 
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Figure 3.7 Spectral radiance of a dark cumulus cloud 


(From Ref.11) 
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Figure 3.8 Spectral radiance for high altitudes (From Ref.6) 
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resembles that of a blackbody with the cloud’s temperature 
except in the region of 5-8 um. [In Figure- 3-10-we have -the 
integral emittance E, transmittance T and absorption A for 
cloud layers with mires sc 6z in meters. 

On a cloudy day we have large variations of the average 
radiance because of the various ways the clouds can scatter 
the sunlight. For a water droplet the scattering is much 
greater in the forward direction than to the side. For that 
reason a thin cloud near the sun seems brighter than it would 


if it was further away. 
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Figure 3.9 Spectral emissivity of cloud layers thickness 62 
(From Ref.11) 
This phenomenon occurs mostly at the short wavelengths with a 


peak at about 3 um. 
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Figure 3.10 Spectral emissivity of cloud layer thickness 56z 


(From Ref.1i1) 


E. MARINE BACKGROUND RADIANCES 
1. Introduction 
The surface of the earth is covered by a great 
amount of water which has some unique optical properties 
Since it is a poor emitter of radiation and a good reflector 
in the optical region of the spectrum. In the infrared 
region, though, water becomes a good emitter under certain 
circumstances. If the water were not disturbed by the waves, 
it would reflect the light like a mirror but, because of the 
wave pattern, which creates various surface elements, we have 
different orientations depending on the slope and height 
distribution of the waves. 
In the IR region the reflection abilities of the sea 
water remain amplified by the good emission of the infrared 
radiation when the waves begin to break [Ref.6). In the 


special case of foamy breaking waves we have a blackbody 
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behavior with a higher radiance for a rough sea than a calm 
one. 
Generally the marine background radiance depends on the 
following parameters [Ref.5]: 
a. Material properties of the sea bottom, 
b. Temperature distribution of the sea’s surface, 
c. Various optical properties of the sea water, and 
d. Wave slope distribution and surface geometry. 
a. t a ° ° otto 
After penetrating into sea water, the solar radiation 
can scatter from the various particles in the sea water and 
from the bottom materials. This scattering occurs only in the 
ioe tO 3.0 tm region falling to the 0.4 to 0.7 um region for 
pure water. This scattering is not important for the eae erie 
of this thesis. 
3. Temperature distribution of the surface 
The temperature distribution of the sea surface, 
ranging from 29° C at the equator to 0° C in the arctic 
regions, plays a major role in determining the total radiance 
of the sea. That radiance can be further effected by the 
presence of currents such as the Gulf Stream, giving 
variations of several degrees Celsius, but mostly is 
determined by the temperature variation from point to point. 
From Figure Soe Ret... 2 | we see a variation in 
temperature of 0.6° C in meaSurements made under evaporative 


conditions at the surface relative to measurements made one Ii 


o1 


mm below the surface. That slight cooler variation “oi@rene 
0.6°C is due to radiation exchange rate of evaporation of the 
heat flow below and above the surface [{Ref.5]}. 

4. t erties of sea water 

The top layer, 0.01 cm of thickness, of sea water 
most determines the thermal radiance [{Ref.5]. The subsurface 
part of the water does not create any radiation. At 
wavelengths above 3 um sea water is opaque. At wavelengths 
from 2 to 4 um the thermal emission of the surface is less 
than the reflected radiation which dominates the wavelength 
range covering 99.02% of the sun’s total radiant energy. 

That leaves only 0.08% for the 8 to 14 um region which is 
mostly effected by the thermal radiation. In Figures 3.12 
through 3.14 (Ref.5] we have the emissivity and reflecuee 
vs angle of incidence calculated from average data from 2 to 
15 um, the reflection from sea water vs. wavelength and the 
incidences of refraction vs. wavelength of the sea water. 

5. Wav t ut ur t 

The surface geometry and wave slope distribution as 
presented by Cox, and Munk in Ref.12 play a major role in the 
refraction and reflection of the incident radiation on the 
sea surface. 

In Figure 3.15 [Ref .5] we have the reflection of the 
solar radiation from a flat surface with o=0 and from a rough 
surface, caused by a wind of 4 Beaufort (11 to 16 knots), 


With o=-0.2 (where of represents the square root of the 


OZ 





variance of the wave slope distribution). From the figure we 
can observe that the albedo R on the vertical scale varies 
starting from a value of 0.02 when the sun is at the zenith 
(¢=0°) up to 1 when the sun falls at the horizon [g=307 }e for 
a flat sea. For the rough surface multiple reflections and 
shadowing begin to be effective causing the reflectance to 
fall to 20% near the horizon and the emissivity to rise to a 


value of 80%. 
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Figure 3.11 Thermal structure of the sea boundary layer 


(From Ref. 5) 
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Figure 3.12 Reflection of sea water surface at 0°,60° and 


80° angle of incidence (From Ref.5) 
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Figure 3.13 Reflectance and emissivity of water vs incidence 


(From Ref.5) 
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Figure 3.14 Index of refraction of sea water calculated 


from refractivity data (From Ref.5) 
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pauore 3.15 Reflection of solar radiation from a flat and 


from a rough sea surface (From Ref.5) 
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IV. THE BACKGROUND MEASUREMENTS AND ANALYSIS PROGRAM 


A. GENERAL 
1. Introduction 
The BMAP code is a program being developed by the 
ONTAR corporation as a supplement to the NSWC clutter code 
for use with the Navy’s fleet defense systems under the 
requirement for detecting the infrared signature of a target 
against a clutter background. To achieve this requirement 
nigh spatial resolution data will be needed from various 
infrared sensors on airborne or surface based niattoume over 
a wide area which will provide the required threat warning, 
acquisition and target tracking. The total field of view 
(TFOY) ase oe searched continuously and, in combination 
with the small pixel size instantaneous field of view, 
results in a large throughput for the signal processors is 
needed. This will cause false alarms originating from the sea 
clutter and clouds that will put some limitations on the 
performance of the system. 
For this reason in support of these defense systems the 
background measurements and analysis program was developed in 
support of the NSWC clubtier ceae with the following 


capabilities [Ref.13 and 14]: 


a. Analysis and modeling of various background scenes. 
py Acquisition of high quality background data 
cn Characterization of the IR scene’s background clutter 
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The program model in including the NSWC clutter code is used 


in the following areas [Ref.14]: 


a. 


Collecting data: The program is to operate under 
field site conditions in a small IBM portable 
computer using real time data for the appropriate 
real time scenes and displaying that data on the 
screen as an infrared scene. 

Data analysis: The option of data analysis, and 
characterization of the data and other cloud scenes 
is available with respect co the atmospheric 
conditions and infrared real time scenarios. 
Modeling: The modeling option of the data and the sea 
clutter is also available, within signal processing, 
ror cabewlatineg the radiance using Monte Carlo 
methods, irradiance and transmittance using LOWTRAN6 
modeling 


Model Parameters (Present Version of BMAP) 
For the current edition of the BMAP [Ref.16] the 


following assumptions have been’ made: 


ale 


Response of the instrument: In the present version of 
the BMAP .the bandpass response for the BMAP 
instrument function in modeling is determined by the 
user. Future modified version of the code will 
determine the instrument response according to the 


wavelength. 


oa 


Geometry of the Clouds: The current edition of the 
BMAP assumes that the clouds under consideration are 
semi-infinite and variable cloud thickness can be 
determined using the Monte Carlo statistical method. 
Wavelength Bandwith of Operation: The present source 
code is limited to the wavelength bandwith region of 
3 to 5 micrometers. The future edition of the BMAP 
will extend that region to the 8 to 12 micrometers 
wavelength band. 

Scattering Model: This model calculates the angular 
Gist ri bution ) fom the clouds by using the phase 
function presented by Henyey and Greenstein [{Ref.20] 
and various other parameters such as the Spectral 
Albedo and the Scattering Asymmetry. Zachor and 
Shette ([Ref.21] determine some of these parameters 
for an average rain radius of 4 micrometers for the 
water droplet distribution used in the present 
version of the BMAP. This study [{Ref.21] can be used 
for various droplet distributions, ranging from Zeto 
32 micrometers. Some of these results will be 
implemented in the future modified edition of the 
BMAP. 

Geometrical Aspects of View: This determines whether 
the cloud is side lit and the scattered radiation 
creates the illuminated surface, or whether the cloud 


is backlit and the radiation from the sun, which is 
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behind the cloud, is passing through the mass of the 
cloud. The present version of the BMAP program 
assumes that the observation is being made from the 
ground; however higher altitude Seeiepansame are also 
permitted. 
£ . Irradiance of the Sure The source code uses the 
transmission coefficient (determined by running the 
LOWTRAN6 program) multiplied | by the irradiance 
existing at the top of the atmosphere to yield the 
radiance at the cloud. 
3. Computer Model for the Cloud Radiance 
This is a program for the 3 to 5 micrometers 
wavelength region which calculates the cloud signature and 
makes estimations for the expected signal - Pees value 
feet. 16). Other calculated factors are the atmospheric 
transmission coefficient using the LOWTRAN6 program, the 
solar radiation scattered by the various cloud and 
atmospheric scenes and finally the spectral response of the 
instruments. Further explicit calculations include the 
scattering of the sun rays according to various geometric 
aspect views (top, bottom, left and right views). For this 
purpose the Monte Carlo statistical method [{Ref.18] has been 
implemented since it is the only one with the desired 
flexibility to be applied to every cloud scene. im addatiron 
to the Monte Carlo approximations multiple expansions are 


used including the doubling method [{Ref.17] and the two 


a 


Stream diffusion models for radiation transport [Ref.19]} for 
the n-stream approximations. 

The main problem of the previously mentioned 
approximations is with the realistic eee dimensional cloud 
scenes. These problems are partially solved, as described in 
chapter 3 of this thesis, by tracking the scattering of the 
individual photons for various scattering phases with the 
obvious disadvantage of anemGnine too much computer time since 
the number of tries for different photons is high. 

In the more specific case where approximated values for 
the radiance model are needed the number of trial photons can 
be reduced. Further indication of the accuracy of the 
calculations is provided by the program by giving the 
statistical uncertainties according to the ike Ee geometries 
and the scattering model. 

By combining different functions such as the spectral 
scattering albedo, the solar spectral irradiance, the 
scattered Signatures spectral bandpass, the instrument 
bandpass, the spectral emissivity and finally the Lowtran 
transmittance the in-band intensity can be calculated. The 
selection of the wavelength for each photon in the scattering 
process can also be determined by obtaining the spectral 
weighting function after combining the instruments components 


and the solar transmittance. 
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B. CODE DEVELOPMENT 


moe present version of the BMAP is the initial edition 
which is based on the modeling of a structured infrared cloud 
background according to various sc sreriic conditions and 
infrared local scenes. Future editions of the BMAP will be 
based on the development of a cloud structure model which 
will be accessed according to the meteorological conditions 
for the local area of observations. This structural model 
will be in conjunction with the frequency and spatial 
BmeractLerization which will be accessed through statistical 
correlation as shown in Figure 4.1 [Ref.16]. The reason that 
this can be done is that we can parameterize a cloud 
structure with respect to its infrared image in terms of 
various parameters such as spatial frequencies, edgeness, 
contrast, etc. 

This can be done by building a data bank of the above 
mentioned parameters for different meteorological conditions 
at the exact moment of observation of the data and in such a 
way that an immediate correlation can be made between the 
input conditions and the data. For this to be achieved 
several cloud images must be analyzed, stored and correlated. 

Finally all that the user will have to do is to input the 
local conditions of observation which will be combined with 
the spatial characteristics for the infrared scene, and a map 
of the cloud with the correct structural features will be 


produced. 
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Figure 4.1 Structural outline of the BMAP program 


C. LOWTRAN6 PROGRAM 

There are three sub-programs inside the BMAP program 
which can be accessed through the function keys of the main 
code, as will be explained in chapter 5 of this thesis, or 
can be run individually as independent programs from the 
directory. 

The LOWTRAN6 is a Fortran program which is designed to 
calculate radiance and atmospheric transmittance t(j} ) 
averaged over 20 cm-! in steps of 5 em-!. The atmospheric 
transmittance [{Ref.1] can be expressed as: 


tT(A) = Ec rCADR ae 
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Figure 4.2 Various Atmospheric models in use with the 


LOWTRAN6 code (From Ref.4). 
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where r(A) is the extinction coefficient (equal to the sum of 
the scattering coefficient t( X ) and the absorption 
coefficient o(A)) and R is the range. 

For the BMAP code further development of the LOWTRAN code 
include calculation of the transmission of Solar Irradiance 
and of the radiance with or without scattering present. 

In the model we have a _ single parameter band model for 
molecular absorption. The code also includes the continuum 
absorption in the 3 to 5 micrometers band by the water vapor 
and the nitrogen and the absorption in the 8 Youu i 
micrometers band by the water is included in the code. These 
two atmospheric windows in the infrared are caused not only 
by the water but also by other absorbing gases like carbon 
dioxide ne 4.3, 15 micrometers), the ozone (4.8, 3.2, 6.3 
micrometers) and the nitrogen. Nitrogen absorption has its 
greatest value in the 8 to 14 micrometers region. 

In the code the earth’s curvature and the refraction have 
also been taken into consideration. 

In the propagation code the atmosphere consists of 33 
layers with variable width (1 Km up to an altitude of 25 Km, 
5 Km up to an altitude from 25 to 50 Km, from 50 to /7Giiae 
and from 70 Km to 100 Km). The user has the choice of 6 
atmospheric models [Figure 4.2] which are the 1962 US 
Standard model, the Mean Tropical model, the Mid Latitude 
Summer model, the Mid Latitude Winter model, the Sub-Arctic 


Summer model and finally the Sub-Arctic Winter model. 
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Furthermore the user has the choice of six aerosol 
models; these are the Rural model with visibility 5 Km, the 
Reeds With visibility 23 Km, the Maritime model with 
visibility 23 Km, the Navy Maritime model with variable 
visibility, the Urban model with visibility of 5 Km, the 
Tropospheric model with visibility 50 Km and finally a user 
defined model with visibility 23 Km. A major difference from 
the LOWTRANG6 program in the IBM 3033 computer of the Naval 
Postgraduate Schoolis that the Advection fog and the 
Radiation fog models in use are not incorporated in the BMAP 
program. 

For the rain cases the ralnaineucea, (extinetion 
coefficient has been defined by the Mie Scattering theory and 
the Marshal-Palmer raindrop size distribution [Ref .4]. 
Appendix A of this thesis includes a table of the cards used 
by the BMAP LOWTRAN code. Further information about the 
structure of these models and about the mathematical 
definitions governing the transmission inside each of these 
models is not within the scope of this thesis but it can be 


found in detail in References 4 and 15. 
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ve CO FOR T 2 TION 


A. GENERAL 

The BMAP program is a code designed to operate in a small 
IBM portable AT computer using the INTEL 80287 math co- 
processor chip and the Tecmar high resolution graphics card 
with an optional hard disk and a 9 track tape drive. 

The code comes with three BMAP Bedford files named BEDFORD1, 
BEDFORD2, BEDFORD3, each 12160 bytes long, which can be 
accessed from the program using the Fi key. 

~All of the data and the text input are accomplished using 
the computer keyboard (standard type) with certain tasks done 
by the arrow and function keys. 

The computer code has a variety of capabilities which can 
be accessed using either the keyboard or a penpad digital 
tablet or a combination of them. The future user can find 
more detailed information in reference 16 about the use of 
the keyboard or the penpad. 

The source code has three main options which are 
represented by three sets of icons. These sets are the 
DISPLAY module, the.ANALYSIS module and the MODELING module. 
A summarized description of these modules is presented in 
figure 5.1 {Ref.16]. Each module set consists of a variety of 
icons which always appear on the left side of the screen ina 
manner corresponding to the F1...F10 keys of the keyboard. 


These modules have one thing in common regarding the 
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2D Plot 
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Histogran 
Edge2 


Metrics? 
2D PSD 
Filter 


Signal 
Processing 
Spectral3 
Spatial2 
Lowin 
Lowtran 


Lowplt 


CODE GAPASLILCITGIES 


DESERTS TION PRINTER . 
OUTPUT 

Access icon set to display new data 

Puts Cursor on the screen & prints 

row, column & pixel value 
Changes the range of pixels displayed 
Generates an on-screen plot of a BMAP 

channel Xx 
Generate an onscreen 3D surface plot 

in 4 perspectives . X 


Displays a histogram of the BMAP Image 
Computes & displays edee operators on 


BMAP imege . X 
Computes passive image metrics 
Computes 2D PSD of data . X 
High & low sass PSD filter & inversa 


2D FFT to compute a filtered image 


o< 


Access icons {£3¢ Signal processing 
Computes Monte Carlo radiance model 
Model cloud spatial structure 
Generates Lowtran input & graph file 
Lowtran 6 


Outputs Lowtran graph X 


Summary of DISPLAY, ANALYSIS and MODELING modules 


(From Ref.16) 
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function keys. The F1, F3, F5, FiO keys “have the Sane 
function no matter what module is running at the moment and 
are: 
Fl: Gets the DISPLAY module 
F3: Gets the ANALYSIS module 
F5: Gets the MODELING module 
F10: Exits the program and returns the user to DOS 
The functions of the rest of the buttons for each module 


are presented in Appendix C of this thesis. 


B. COMPUTER MODULES 
1. Display Module 
This module appears on the computer screen when the 
pregram is loaded and consists of the infrared image on the 
top of the screen with an empty space below (in which the 
different program options appear with a set of icons having 
the following functions/capabilities: 
a. Icon 1 (Fl key): We have No change since the DISPLAY 
module is already present. 
be Icon 2 (F2 key): Allows the user the selection of a 


new image and returns him to the data directory 


icons. 

ey Tcon 3 (F3 key): Brings the ANALYSIS module on 
screen. 

d. Icon 4 (F4 key): Creates a cursor on the infrared 


image on the top of the screen which will provide 
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data information about the specific point under the 
cursor. The data is presented under the image first 
by the number of the channel which that point 
represented when the data was taken, then the sample 
number and finally the radiance (number of counts). 
Icon 5 (F5 key): This key brings up the MODEL module. 
Icon 6 (F6 key): This key makes the adjustment in the 
range of the count values on the display which is 
usually within some specific range of values. The 
program asks the user for an upper and a lower 
threshold value for the new range which will improve 
the ability to redistribute the shades of gray in the 
new range and to improve the clarity in the picture 
of the image. 

Icon 7 (F7 key): When this option is selected the 
program provides a two dimensional plot of one 
channel only. In the beginning the user must provide 
the code with the channel number (from 1 to 16) that 
he wishes to plot. Then the program will provide the 
counts for the required channel anda plot of the 
sample. Each count signifies radiance of 5.6x10-8 
watts/cem@-steradian. In Figures 5.2 and 5.3 [({Ref.16] 
we have a representation of the channels #2 and #6 of 
the BEDFORDI file. 

Icon 8 (F8 key): This key will give the counts versus 


samples for all of the 16 channels and a three 
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Figure 5.2 Two dimensional plot for channel #2 of BEDFORDi 
(From Ref.16) 
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Figure 5.3 Two dimensional plot for channel #6 of BEDFORD1 
(From Ref.16) 
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Figure 5.4 Three dimensional plot of BEDFORD! file 
(all 16 channels), (From Ref.16) 
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dimensional plot of the image, and provides the user 
with different options regarding the way the plot is 
to be observed (view aspects). These options are 
Peon sro DACK == rroant “Lo left, Lert to right and 
finally back to front. Figure 5.4 [{Ref.16] gives the 
three dimensional plot of the BEDFORD1I file for all 
16 channels. 

Icon 9 (F9 key): This key will give the directory of 
the disk and is Similar to the DOS command 
Porreerory . 

Icon 10 (F10 key): This key exits the program and 
returns to DOS. 

ANALYSIS module 


This is the module that provides modeling of the data 


and various cloud scenes with the option to present and 


filter 


the power spectral density (PSD) of an image. This 


module consists of 10 icons with the following functions: 


a. 


Icon 1 (F1 key): Brings the DISPLAY module and its 
functions on the screen. 

Icon 2 (F2 key): This function will give a histogram 
of an image with the number of counts.) on the bottom 
of the screen. (The vertical direction signifies the 
number of pixels for the specific count value. 

Icon 3 (F3 key): This brings up the ANALYSIS module 


and its functions on the screen. 
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Icon 4 (F4 key): This function is intended for future 
use with the revised edition of the BMAP and will 
compute and display the edge operators for an 
infrared image. 

Icon 5 (F5 key): This brings up the MODEL module and 
its functions on the screen. 

Icon 6 (F6 key): This function is also intended for 
future use in the revised edition of the program and 
will compute passive image metrics. 

Icon 7 (F7 key): This will compute a two dimensional 
power spectral density (2D-PSD) for an infrared scene 
and display it on the screen. This is done by first 
calculating the Fast Fourier Transform of the image, 
Lis chaes has not already been done, and then 
presenting it on the screen in a |logarithmic or 
linear scale according to the user’s wish. In the 
case that the linear scale is chosen as an option a 
further option of clipping the DC value is available. 
In Figure 5.5 [Ref.16] we have a two dimensional plot 
of BEDFORDI file with the vertical axis representing 
the power spectrum (linear scale), the horizontal 
axis the sample counts andthe transverse axis the 
channel number. 

fcon 8 (F8 key): This function will provide filtera ae 
of the power spectral density and then display it on 


the sereen by first computing the Fast Fourier 
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Transform of the image and then filtering according 


to the cutoff index. Four- 4ypes of filtering are 


available: 
ie In the sample direction / High pass filtering 
ee In the sample direction / Low pass filtering 
2). In the channel Girecvion y High pass 
filtering 
4. In the channel direction / Low pass filtering 
8 
3 DAYS 8 HRS S2 ri | ae 
els SEE ial 


LW IPOREAR SPECTRUM) 





BEOFORO DATA 
TAPE? 





Figure 5.5 Two dimensional plot of BEDFORD1 file for power 


spectrum (From Ref.16) 
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When the filtering is occurring in the channel 
direction the cutoff index can be a channel from 1 to 
8. When the filtering is in the sample direction the 
cutoff index is a sample number. 
The roll off coefficient C then must be inserted 
found from the formula [{Ref.16]: 
Roll off = e-D%vc (5-1) 
where D=distance from the cutoff, and 
C=roll-off cocmeicicne- 
This will give the rate with which the signal falls 
offi vafiter the lcuron ie 
The previously mentioned filtering techniques can be 
used individually or all together. Then the program 
will perform the inverse Fast Fourier Transform of 
the image and will display it on the screen. 
in Icon 9 (FY key): This option is intended for future 
use and will access further icons from the code for 
Signal processing. 
Te Icon 10 (F110 key): Exits the program and returns the 


user to DOS. 


3. MODELING module 
This module of the BMAP calculates the radiance using 
the Monte Carlo method, previously described in chapter 3 of 
this thesis, and the reflected and transmitted intensity for 


each cloud layer ina cloudy infrared scene. A further option 
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Se this module includes the LOWTRAN6 code with its 
accompanying subprograms LOWIN and LOWTRAN plotting. More in 
detail the module consists of the following icons: 

a. Icon 1 (F1 key): Brings up the DISPLAY module and 
it’s icons on the screen. 

oy Icon 2 (F2 key): This function will compute the 
radiance model using the Monte Carlo statistical 
method model. That model uses values taken from the 
LOWIN subprogram which can be changed according to 
the desire of the user. That means that the LOWIN and 
LOWTRAN programs must have been run ahead before the 
radiance is calculated. Otherwise the values for the 
Monte Carlo statistical method analysis will be from 
a previous run and the results will be wrong. The 
spectral modeling icon then will calculate the 
reflected and transmitted intensity for each cloud 
layer. It must be noted though that this function is 
yet not fully implemented. 

Cc. Icon 3 (F3 key): This brings the ANALYSIS module and 
ivsectunecrt ions On Ssereen. 

d. Icon 4 (F4 key): This is intended for future use in 
the revised edition of the BMAP and will give a cloud 
model Spatial Structure. 

e. Icon 5 (F5 key): Brings the MODEL module and its 


functions on the screen. 
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leon 6 (F6 key): This function is the laa 
subprogram for the LOWTRAN6 program. Its use is to 
insert the values required to run the LOWTRAN program 
and the LOWPLOT program. This subprogram consists of 
4 cards, presented in Appendix A. Further explanation 
of the function of these cards and the various. 
variables that can be inserted can be found in the 
Chapter 4 of this thesis and in Ref.15. 

Icon 7 (F7 key): This option will use the values 
which are assigned to the Lowin subprogram (icon 6) 
and run the LOWTRAN program. More details about the 
LOWTRAN6 program can be found in Chapter 4 of vais 
thesis and in Ref 15. 

Icon 8 (F8 key): This option will give the output 
LOWTRAN plots. When this function comes to the screen 
it will present a card, an example of which is 
presented in Appendix A, with different variables 
according to the type of lines that are going to be 
used in the plots. More details can be found in 
Ref .1i5. 

Icon 9 (F9 Rees Wate will bring up the directory 
from the disk and it has the same function as the 
“Directory” command of DOS. 

Icon 10 (F10 key): This will exit the code and return 


the user to DOS. 
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As a general remark for all modules it should be noted that a 
help routine is always available to the user which can be 
accessed by pressing the <Alt> key in conjunction with one of 
the Fi...F1i10 keys. Tho code then will give a short 


description in the function of the specific key. 
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Vix BMAP TAPE FORMAT 


A. GENERAL 

The BMAP code consists of three data files called 
BEDFORD1I, BEDFORD2, BEDFORD3 with the name originating from 
the place of data measurement on top of a hill at an 
elevation of 250 feet above mean sea level in Bedford, MA 
from 10 to 14 September 1984. 

The experiments were made with the participation of NSWC, 
NADC, NRL and the Raytheon Corporation in a hole in a solid 
cloud deck in the sky that was backlit by the sun. The sensor 
was a Raytheon dual band scanning radiometer in the 
wavelength spectrum from 4 to 5 and from 8 to 11 micrometers. 
The sensor had two 16 element arrays with an instantaneous 
field of view (IFOV) 0.33x0.33 millirad and Total Piet. 
View (TFOV) 16x512 pixels. It must be noted that the arrays 
were boresighted and calibrated with reference to a two 


temperature internal blackbody standard. 


B. INTERNAL CONFIGURATION 
1. Introduction. 
The BEDFORD data files are arranged according to an 
interactive computer code called ATITP based on the NRL code 
6520 implementation of the NATO standard Image transfer 


[Ref.24] requirements. The computer code gives as an output a 
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9 track magnetic tape which can be used for infrared data 
ianerchange. 

The NRL code 6520 is a VAX/VMS 11 disk structure which 
will be imposed with an internal file structure in such a way 
that the information can be transferred to a 9 track magnetic 
@ieemtor further data manipulation. Its purpose is to enable 
the various users to access the data inside the NRL magnetic 
tapes. 

The BEDFORD files in the BMAP are arranged in the same 
way in a format consisting of three parts as in Figure 6.1 
[Ref .24]. 

The first part has the name “Header 1" and is a record of 
the source of the file data BEDFORD files anda brief 
description of the data in part two and three. "Header 1" 
also gives information about the origin of the files, the 
Gata organization and the structure of the files ina 
standard format. 

“Header i" consists of three basic parts. The first part 
consists of Pheer tes irom 1 aoOmce Genta yeaeerae Origin of 
the image. The second part consists of the bytes from 25 to 
80 giving a description of the rest of the data in the file. 
The third part uses the bytes from 81 to 104 giving 
information about the x and y sizes and the type of the data. 
Finally we have the bytes from 105 to 128 which are 


unassigned at present and will be used for expansion and 
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modification for future versions of the code by specifying 
multi-image -fites.-- 

The second part has the name "Header 2” andis a tape 
record of the structural detail of the data following in part 
three. 

With "Header 2” we can have free format documentation of 
all the necessary parameters for the data such as the 
location of the significant image features and documentation 
relating the recording environment. 

It must be noted that at the end of each header there 
follows an "End of File” (EOF) mark which usually creates a 
problem since it makes each header a separate file, 
consequentially making the main file (BEDFORD file) 
unreadable without some special decoding technique. The 
problem can be eliminated by ignoring the headers’) and 
skipping the first two records providing, of course, that we 


know the image data structure that follows. 


-|Header 1 Fixed 128 Bytes 
[Header 2 







Header 2 Fixed 128 Bytes 


Data Record 1 Fixed in Header 1 





Data Record 2 





Data Record N 


Figure 6.1 BEDFORD files byte structure (From Ref.16) 






Fixed 128 Bytes 
Fixed 128 Bytes 
Fixed in Header 1 
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The third part is the data of the infrared image, 
described in this chapter introduction, stored-as one channel 
of data per record. 

The data file representing each image in the tape record 
differs in form from the organization of pixel data recorded 
from the measurements sensor [Ref.25]. Thus we can have each 
image written on the tape as a separate file without 
necessarily having each record containing one image line. The 
form of this configuration permits the user to use the format 
to transfer very long image lines (16 channels for the 
BEDFORD files images) using a small tape buffer size of 4K 


existing in small computers, such as the IBM AT, and in mini 


computers, such as the Masscomp computer of the IRSTD. 


UISIAT INIRILI 
8i5t | 13) [2/3) 
Tape ee le 







ORIGINATOR CODE 

DATE TAPE MADE 

NRL INTERNAL FILE IDENTIFIER 
NUMBER OF RECORDS IN FILE 
SIZE OF FILE 2 

TWO BYTE INTEGER URLUES 

ONE ENTRY PER SAMPLE 

TWO BYTES PER “CHANNEL” 
RERL MANTISA BYTES (UNUSED) 
REAL ERPONENT BYTES (UNUSED) 
SAMPLES PER LINE 

LINES PER IMAGE 

NONCOMPLEH INTEGER DATA 

16 UNUSED BYTES. 






Figure 6.2 BEDFORD file “Header 1" configuration 


(From Ref.16) 
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The basic format of the data is based on the tape 
Character (byte). The same format can be used for 7 track 
tape as 6 bits and in 9 track tape, as for the BEDFORD files, 
at 8 bitse 

a: t of “#H e 

The "Header 1 field consists, of 128 byteSmuamers 
eight byte fields in a configuration according to Figure 6.2 
[Ref .24]. 

“Header 1° for a 7 track tape is coded in BCD code even 

parity and for a 9 track tape in ASCII code odd parity. 


In the “Header 1" we can find the following information: 


Bytes Presenting Information 
too Identity of the originator with 4 symbols for 


the country and 4 symbols for the laboratory 


J vomrG Date of data: Year, Month, Day 

Li tOw es Name or number of file, image, or program 
Unies 

22 E@ see Number of records per file (R) 


(R=N+2 to account for 2 Header records) 

33° temac Number of tape characters in "Header 2" (L2) 
with the requirement L2>0 

al to 74s Number of tape characters per image data 


record (L3) or Number of bytes in the program 


units (L3=80n With = ieee 
49 toGgoe Number of entries per sample (C) or Number of 
program units in the file (one channel 
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complex data are treated as two channel data 
with C=2 x number of channels) 

57 to 64 Number of tape characters per channel for 
integer data (I) or Number of card images on 
the file (I=0 for real values) 

eomto f2 Number of tape characters for mantissa of a 
reali valued channel (Em) (Em=0 for integer 
values) 

foe to 80 Number of tape characters for exponent of a 


real valued channel (Ee) (Ee=0 for integer 


values) 
eo LO 88 Number of samples per line (S) 
Be to 96 Number of lines per image (L1) 
97 to 104 Type of data values (T) with: 


T=0 for non-complex, integer 

T=1 for complex, integer 

T=2 for non-complex, real 

T=3 for complex, real 

a4 for source code card images (80 

characters per source code line) 
mea to 128 Currently not assigned 

3. Structure of “Header 2” 
The file “Header 2" in its structure defines a series 

of fields containing ASCII character strings taking advantage 
of the lack of specification of the NATO standard format. In 


those fields the numeric string is stored and the binary 
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values are represented by the letters "T"” (true) or "F" 
(false) according to the Fortran standard. The "Header 2” 
structure format is shown in Figure 6.3 [Ref.24] with each 
square representing one byte. 

The “Header 2" configuration contains a description of 
the image including the physical significance of the various 
channels, meteorological conditions of the day when the data 
were taken, description of the scanner used, location of 
objects in image, etc. Its length is variable between 0 and 


4096 characters [Ref.25] and has the same coding as the 


“Header 1”. 
olo} .42 A4 HEADER 2 VEASION 
x XIXIXI XIX IX A352 ORIGINAL DATA TIME 
X{ x] xixix{x1xlx COBE 
a X/X/X/X/x]x 
X|x{ x] Xi x/x1/x1x 
| | Li TRUE IF LONG WAVE DATA 
ie L2 TRUE IF MID WAVE DATA 
Pf {2 CORRECTION PGM. CODE 
Se 14 CORRECTION P&M. VERSION 
cer re 14 DAH DATA FORMAT 
ae ee 14 FIRST CHANNEL OUTPUT 
ae ccee 14 LAST CHANNEL OUTPUT 
Pee 14 FIRST EDEN DATA POINT 
SESS | 14 LAST EDEN DATA POINT 
ry toy 14 FIRST ODD DATA POINT 
Bees. 14 LAST ODD DATA POINT 
She “i 64 BYTES FOR FUTURE 


EHPANSION 


Figure 6.3 BEDFORD “Header 2” configuration (From Ref.16) 
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In the file the first 4 bytes are used to specify the 
header version number which can always be changed according 
to whether new information has been added or removed. At 
present the version is 00.2. The next field is 32 bytes long 
and is used to specify the time at which the data were taken. 
We then have two 1 byte logical fields that signify the wave 
band of the data. [f the letter "T" is on the third field 
then the data are longwave (8-12 micrometers). If on the 
other hand the letter ‘"T" is on the fourth field then the 
data are midwave (3-5 micrometers). For the BEDFORD files the 
letter “T" is on the fourth field. 

The next bytes are used for file correction. According to 
NRL standards the usual method of doing the correction of a 
mle iis by using two data correction programs called 
“CORRECTIN" and "ERRFND" £a[Ref.24]. The program used is 
symbolized in the "Header 2" file in the next field by the 
two digit code "1" if the correction program is the 
Beonnme@liN and the two digit code 2 if the correction 
program is the "ERRFND". For better insurance the version 
number of the correction program is stored in the next four 
bytes position. 

The next field, four bytes long, is used to give the VAX 
tape data version number from which the data were taken. 

The next six fields, also four bytes long, are used for data 
field editing information as an insurance for a trace back to 


the original form of the data and contain the last channel 
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output, first even data point, last even data point, first 
odd data point, and last odd data point. Finally the last 64 
bytes are used for future expansion. 

4. Data structure 

The data inside the BEDFORD files are stored in fixed 
length records the number of which can be found from 
information given by the “Header 1". For the BEDFORD files, 
Figure 6.2, the number of records is 18. In each record we 
have one scan line of data. That means that for the 18 
records we have 18 scan lines of data for each BEDFORD file. 
The image data are binary odd parity and one’s complement 
[Ref.25] in integer format which remains the same for all of 
the other types of format as described in “Header a 
configuration for the bytes 97 to 104. In the integer fommat 
the most important part is recorded as a first tape character 
with the most significant bit (the farthest left) being the 
Sign bit (QO for positive and 1 for negative). 

For the real format the exponent and the mantissa, which 
are base 2, are reduced to integers and are then stored as 
separate integers (for each sample) on the tape. For the 
complex format the integer values (T=1 in bytes 97 to 104 of 
“Header 1") are stored as integers and the real values (T=3) 
are stored as real values. We then have two values for each 
channel [{Ref.25]. The first one is the real part and the 
second one is the imaginary part as in Figure 6.4 [{Ref.25]}. 


Each channel in return consists of a real and an imaginary 
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part and each sample consists of 2 channels (2 real and 2 
imaginary parts). 

The NRL data banks have various infrared images in NATO 
standard format in various computer formats like the 
following cases. 

In the case of multi-channel data the register and the 
non-register parts are stored as separate images in separate 
files. 

Another method for the multi-channel case is to store the 
data on a pixel interleaved basis and to remove the values of 
the different channels as adjacent bytes for one sample as 
in Figure 6.5. This figure shows storing methods for three 


color non-complex data and two channel complex data. 


2 channel, complex data 


Channel i Channel 2 Channel 3 
a : be | . j : i‘. 





Sample M Sample Ml 


Figure 6.4 Complex Format for Image data 


Each sample (M, Mi) consists of 3 channels. For the case 


of the BEDFORD files the data are written according to the 


standard NATO notation as binary and one’s complement which 
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is the same as the two’s complement since all data are of 
reasonable magnitude and positive. 
In Figure 6.6 [Ref.24] we have a representation for the 


data format structure of the BEDFORD files. 


3, Colonsonon=conmpillexwscaca 


ax | ow: cae | | oes Chae _ 6 
Sample M sample Mi 


Figure 6.5 Multi-channel data storing methods 


At the end of the tape a marker <END OF FILE> is placed 
Signifying the end of the data records. After that a new tape 
with data records can be added. The final end of the tape is 
marked with two markers <END OF FILE>. 

5. Physical Interpretation of the data 

The BEDFORD files data were taken, as already 
described in the previous chapter, using a sensor consisting 
of two bore-sighted telescopes measuring Intensity. In front 
of each telescope there was an optical bandpass filter 
covering the sight of sixteen infrared detectors. These 
Infrared detectors were HgCdTe for the longwave band and InSb 
for the midwave band. The first telescope produced i/f noise 


which required data collection Circuits with quasi-de 
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response with a 0.5-1000 Hz and 3 db passband. The electronic 
bandpass of the second one was 0-1000 Hz at 3db [Ref.26]. 

Each telescope scanned across 2.2° of azimuth with a 
speed of 36.0° per second and for a scan time of 60 seconds. 
The intensity then was sampled every 0.096 mR which created 
300 data points per scan. It must be noted that the sampling 
time was less than the dwell time and since the angular 
resolution was 0.33 mR that creates 3.44 intervals of 0.096 
mR per dwell. 

For each detector the random noise level was 2.5 mV. For 
that reason the data digitization was from 0 to 11 (12 bits) 
(00-4096 decimal) so that each count was also 2.5 mV. 

The random noise error was 2-3% since the maximum relative 


intensity variations spanned about 50 counts. 





First Channel Data 
Second Channel Data 






Number of 






® 
» 
e 






| Records 
Third Channel Data Deduced 
e @ e From 
Header 1. 
Last Channel Data Typically 
There are 





18 Records. 


Number of Bytes In a Data Record 
can be found in Header ft. 
Figure 6.6 Data structure format of BEDFORD files 





(From Ref.16). 


a 


VII. AGA DATA FORMAT 


A. GENERAL 

The AGA Thermovision 780 is a thermographic device with 
digital recording and data analysis capabilities able to 
detect infrared radiation in the middle infrared spectrum (3- 
5 micrometers) and/or the far infrared spectrum (3-14 
micrometers) using an InSb detector for the first and a 
HgCdTe for the second. Further information about the AGA 
system can be found in the AGA manual in Reference 28. 

The AGA is connected t0O a micro-computer of West German 
origin, BMC IF 800 model 20, which uses aZ-80, 4 MHsa clock 
microprocessor and a CPM operating system. 

The computer operates with a software package consisting 
of three major programs called DISCO, IMAGE PROCESSING and 
UTILITY which interface with the AGA Thermovision 780 and 
record infrared Images in a 16 bit word form on a diskette in 


one of the systems drives [Ref.28]}. 


B. IMAGE/DATA TRANSFER 

The infrared image in data form can be transferred at any 
point from the diskette into a RAM inside the computer for 
further processing. 

The user though must write a program in either FORTRAN or 


BASIC (provided that they are installed into the QEbs 
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computer) and call a subroutine from the DISCO diskette 
called REDA. 

The subroutine REDA is written in assembly language and 
its purpose is to transfer an infrared image from the data 
diskette to RAM memory inside the computer as an 8192 or 
16384 pixel file, depending whether the form of the image is 
FIELD or FRAME [Ref.27]. Each image inside the RAM will 
occupy 8192 or 16384 bytes accordingly with a 77 16-bit word 
directory data. The number of images that can be transferred 
depends exclusively on the size of the RAM. 

The image is then stored as a matrix of 128 elements 
(width) by 64/128 elements (height) in a sequential method 
starting from the left hand corner of the image and passing 
from the end of the first line restarting at the left most 
image point at the second line of the image. 

The REDA subroutine can be called by writing inside the 
user’s program: 

“call REDA (Paraml, Param2, Param3)" 
where Paraml is an integer from 1 to 38 corresponding to the 
image number, Param2 is also an integer corresponding to the 
address of the first element of the image directory buffer 
and Param3 is a single byte integer corresponding to the 
first pixel of the image buffer. 

The stored image in the RAM has ae specific directory 
stating the conditions and the parameters of the system in 


which the image was recorded. The format in which the 
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disk ie I ascii I high and low byte holds each 
label l I 1 | ascit 
im.no. ." Fall 2 I 16 bit int I image number 
St I 
RANGE i 3 I 16 bit int I 
LEVEL 1 4 I 16 bit int I 
I I I 
abs(A) Et. SG I 24 bit int 1 abs(A)=65536*DI(5)+256*DI(6)+DI(7) 
B: I 8 I 16 bit int I 
abs(C* 1000) I 9,10,11 1 24 bit int I abs(C)=(65536* DI(9)+256*DI(19)+ 
I I l +Di(11))/1000 
_ Deal* 10 1 12 1 16 bit int | calibration distance (m) 
| | dcal=DI(9)/10 
; sign (A,C) Pyle I 16 bit int I bit O: sign of C 
|  aonieasten, oie 
| I I I 
~ EPS*100 I 14 I 16 bit int I emissivity obj. * 100 
DIS 1 15 1 16 bit int 1 object distance (m) 
Tamb* 10 1 16 I 16 bit int [I ambient temperature (K) * 10 
_ Tatm* 0 ane? 1 16 bit int | atmospheric temperature (K) * 1% 
| TRMS* 190 I 18 I 16 bit int I transmission * 100 
' alpha* 100000 1 19,20,21 1 24 bit int I alpha=(65536*DI( 19 )+256* DI( 29)+ 
I I l +DI(21))/190000 
| beta* 100000 I 22,23,24 1 24 bit int I beta=(65536*DI(22)+256* DI( 23) + 
; 7 I I +D1(24))/100000 
I l I 
D/R sw I 25 | 16 bit int I direct/relative switch 
I I I lsrelative ,O:direct 
EPSR* 100 I 26 1] 16 bit int I emissivity ref., epsr=DI(14)/100 
XR } S27 I 16 bit int I horizontal position ref. 
YR I 28 I 16 bit int I vertical position ref. 
VALR ice [ 16 bit int I digital value of reference point 
Tref*10 I 30 1 16 bit int I reference temperature (K) * 10 
I afl I 
SCTY I 3i I 16 bit int I scanner type 
SCVER 1 8232.33, 35 "1 ascii I scanner version 
SERNO Ise I 16 bit int I serial no. 
FILCOD i 936, 37,56 ascil | filter code . 
. FQOVE2 I 39 I 16 bit int I fov=DI(39)/2 
APERT*I0 1 40 | 16 bit int 1 aperture, ap=DI(40)/10 


Figure 7.1 REDA subroutine directory config. (From Ref .28) 
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directory of each image will be stored inside the RAM 
from the REDA subroutine is presented in Figure 7.1. 

With the use of the REDA subroutine the user can connect 
the data generated by the AGA Thermovision 780 with other 
processing software like the IBM operating BMAP program. 

The main problem is the incompatibility of the 16 bit 
image of the AGA to the 8 bit image of the BMAP program and 
the difference in the data format. Further discussion of the 


differences is presented in Chapter X of this thesis. 


Vit. B SSING IN T IRSTD SYST 


A. GENERAL 

The NPS-IRSTD data processing system is based on the 
initial version of the system the IRST (Infrared Search and 
Track System). 

The IRST is a system with automatic detection capability 
for infrared targets oxeeiaea to the system from the 
Electronic Processing Equipment (EPE) after the necessary 
clutter rejection (which is the lower blue sky target signal 
insertion loss with a corresponding high clutter rejection), 
in conjunction with synthetic display capability Gnas 
time basis. 

The EPE consists of the following components [Ref.22]: 

a. Signal Processing Subsystem (SPS) with sub-components: 
1. Background Normalized Unit (BNU) 
2. Signal Buffer and Power Supply Unit (BPU) 
3. Data Conditioner Unit (DCU) 

b. Data Processing Subsystem (DPS) with sub-components: 
1. Running software 
2. DPS hardware 


c. Display and Control Subsystem CoCo). 


B. INTERNAL CONFIGURATION OF THE IRST/IRSTD SYSTEM 
13) §S al cessing 9 m 
The Electronic Processing Equipment Unit presented in 


8.1 and 8.2 with its subordinate units, a full deseri pticnwes 
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Figure 8.1 Signal Processing Subsystem Layout 1 
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Figure 8.2 Signal Processing Subsystem layout 2 


(From Ref .22) 
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which is given ip the IRST manual in Ref.22, receives a pre- 
processed raw video from the optics, spectral filters, cold 
stop, detectors, pre-amplifiers and bandpass filters through 
180 detector channels from the scanning infrared optical 
unit. These channels come from two 90-detector linear arrays 
(in the image plane) separated by ai/2 degree, which are 
spectrally filtered in different wavelength bands. The 
existence of the two arrays ensures the coverage of the 
specified elevation angle with a total field of view of 360° 
maz a mMiith . 

These detector channels (180 total) providing the raw 
video coming from the detector’s pre-amplifiers and bandpass 
filters interface with the Signal Processing Subsystem and 
are filtered through a technique which keeps the specified 
per scan constant false alarm rate (CFAR) with a 
ecrresponding low target insertion loss. The normalized 
detector channel data are then synchronized by the SPS unit 
with the corresponding coordinate data so the Data Processing 
subsystem data are buffered, digitized and formatted. 

Zs or Su tem 

The unprocessed video after the Signal Processing 
Subsystem goes to the Data Processor Subsystem (DPS) where a 
further correlation/track and target discrimination is been 
performed so the IRST false designation requirement is 
maintained. Further on the DPS unit keeps the tracking of all 


the pre-designated targets and reports the track status on 


oo 


the Display and Control Subsystem in the form of target 
amplitude, target coordinates and detection array number. An 
additional ability of the DPS unit includes reporting and 
target designation through the NTDS system of the home 
platioum# 
oe d nt u 
This subsystem provides the necessary control of the 
IRST system by displaying both the unprocessed and the 
processed video in conjunction with the system status. 
Further on the system provides identification of any track as 
anoise or clutter track, target designations, blanking zone 
capability, scan false alarm rate control and other track 
store information according to the user’s request. More 
information about the DCS unit can be found in Ref. 22. 
4. Data Processing 
The data processing functions are performed by the 
AN/UYK computer of the IRST system and by the Masscomp 
computer of the IRSTD system with its software with the 


following modes of operation [{Ref.22]: 


a. Data Input (for reporting detections) 
De External Function (for control of the computer) 
Cc. Interrupt Input (for the time code to the computer, 


ships heading, and reporting of the system status) 
ae Data output (for the testing of the Data Conditioner 


Unit  funeta ons). 
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C. DATA CONFIGURATION 
1. Introduction 

This is the mode in which the signal processor target 
detections will be input. The input data consists of two 
words, each with a length of 16 bits for the IRST system and 
@eera bits for the [RSTD system, which will be transferred in 
parallel for each message. For the 16 bits case of the IRST 
system each input consists of 8 bits for the amplitude, 1 bit 
for the array number, 7 bits for the channel number 
(elevation information), 14 bits for the azimuth degrees and 
the rest are control bits [Ref.22]. The memory words will be 


transferred to the Data Conditioner Unit with a maximum rate 
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Figure 8.3 Format of the 16 bits words input to the DCU 


The threshold settings though will be less than 1K 
words/second with a maximum peak rate of 50K words/second. In 
Figure 8.3 [Ref.22] we have the information format of the two 


Me O1t words. 
nO 


Where: A. Array designation with O-leading and 
l=trailing 
ELEVATION. Detector number from 0 to 89 
AMPLITUDE. Detection Amplitude Magnitude 
AZIMUTH. Azimuth counter value at the time of the 
report to the Data Processor 
E/E. Empty/Full indicator with l=Empty and 02800 

The input queue will be 1024 words in a circular buffer 
which will permit a maximum of 512 two-word unprocessed 
reports to be queued before an overflow happens. The 16 bit 
format though in the IRSTD version will be used for testing 
purposes only. The actual format will be in the 12 bits word 
versions where the 1 to 8 bits will be used for the Analog to 
Digital Convertors (ADC) of the IRSTD (12 ADC’s x 15 max.=180 
channels) and the remaining 4 bits in each word will be used 
for position, frame and rotation synchronization. 

2. Synchronization and Geom of the Data 

The IRSTD data coming from the optical sensor unit 

can be considered to be an array of 180 rows (representing 
180 channels) by 60,000 columns (representing 360 degrees 
SSMS EIS of the optical sensor). 

The data coming from the sensor will then go to the 
Masscomp computer ina form of stream of data words, each 
word representing a data value for each of the 180 sensors in 


the IRSTD scanner unit. That data stream will then be 


a 


repeated 60,000 times for every 180 data values representing 
a revolution of the optical unit. 

The Data acquisition will then have to face three 
problems. The first one will be the acquisition of the 
continuous stream of the data for the time the data stream is 
flowing without any loss of data and consequently of the 
information carried. The second problem will be the 
acquisition of the data for the sort time in which the system 
is supposed to receive the data following with a period where 


the flow is stopped. 


meso) ACO TJ. 2... A(O0,60000 ) 
PO) ACL Ie 6. kos A(1,60000) 






eo 
rows 
Mi@iate)) ACLGO,1)°-. . A(i80,60,000) 
Figure 8.4 Data .matrix representing 180 sensors otf the 


Optmeal Un1t 


The third problem will be the format of the data under 
the assumption that they form an array of 180 rows by 60,000 


columns representing the 360° rotation as in Figure 8.4 and 
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more specifically the determination of where the boundaries 
of the array will be. 

The first problem will be solved by controlling the data, 
forcing them to have the correct rate and input timing by 
synchronization of the data stream. This synchronization will 
be made by two leading bits in the beginning of each word. In 
such a way the software of the Masscomp computer will be 
synchronized with the optical scanning unit. 

The first bit will be set to the “one” state at the 
beginning of the rotation of the optical unit and then to 
“zero state for the period of time that the rotataleue:s 
being made. 

The second bit will be set to one for the first 180 
sensors and to zero for the next 179 sensors. With the nedp 
of the first leading bit the Masscomp computer will be first 
synchronized with the optical unit and then will start 
reading the data until the second bit changes value, meaning 
Lhave = tine first 180 sensor values have been read and 
transferred to the Masscomp computer. 

The second problem will be solved by buffering the data 
stream in two separate buffers. In that way when one buffer 
is filling with the incoming data, the next buffer will be 
available for storing and processing of the data into the 
magnetic tape drive of the IRSTD processing system. The data 
will be parallel input into the buffer every synchronized 


clock pulse as presented in Figure 8.5. 
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The storage in the high speed tape recorder will be in 
the form of blocks according to the time required for that 
block to be stored in the magnetic tape. A special provision 
for the storage will be that the time required to read the 
next block is longer than the storage time so as to eliminate 


the possibility of overflow. 
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Figure 8.5 Transfer of data in the Masscomp computer 


Finally the third problem will be solved by recognizing 
which is the first element (sensor) of the optical unit, 
representing a data word in the 180x60,000 array in Figure 
38.4 by A(0,0). at the top left corner, being at the same time 
the first element of the matrix. That first element is the 


most crucial for the data synchronization later on. 
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TX: DATA ANALYSTS 


A. GENERAL 

In this chapter we will discuss the way the data from the 
BEDFORDI to BEDFORD3 files are decoded from their original 
form to alpha-numerical form and then how the data of which 
they consist were transferred to three output files called 
BEDFORD1.OUT, BEDFORD2.OUT, and BEDFORD3.OUT. Finally we will 
discuss how these output files were transferred to the 
Masscomp computer of the NPS-IRSTD system. 

We will also discuss the method in which the decoding 
computer code was developed explaining each component of the 


code in detail. 


B. COMPUTER CODE ANALYSIS 

1. Introduction 

As it was already explained in Chapter V of this 

thesis the BEDFORD files consist of two header form of 
information followed by a matrix of integer data. Each header 
consists of two 128 byte alpha-numerical (character) strings 
with special and specific details about the data that follow 
in an integer matrix of infrared information in integerx2 
form: 

The whole infrared image, represented by the data matrix, 
consisted of 16 lines of 371 data points making a total of 


5936 data points. 
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At the beginning of the code analysis an attempt was made 
to read the data by conventional methods by examining each 
file with the DOS command "BROWSE". 

The attempt was unsuccessful due to the fact that the 
data in the BEDFORD files were encoded. A computer program in 
an advanced language was required that could read that data. 
We then contacted the ONTAR company, which developed the BMAP 
program to explain the coded form of the data and a possible 
way to decode the data. The ONTAR company responded by 
sending some information about the form of the headers and a 
computer code written in Microsoft Fortran which they claimed 
could read the data. The ONTAR computer code is presented in 
Appendix D. 

We then tried to use the ONTAR code to read the BEDFORD 
files but with unsuccessful results since the Microsoft 
Linker running the ONTAR computer code could not recognize 
mae tunetions  OPENFH”, "“RDEMEM’, “IADDRS" and "CLOSFH" found 
in the code. 

The ONTAR company was contacted once more and asked for 
the above mentioned functions. They responded to us Saying 
that these functions are proprietary to the company and they 
could not be given to us. This response made the code 
provided by them unusable. Another method had to be found to 
read the data from the BEDFORD files. That method was to 
develop our own computer code in RM/FORTRAN that would be 


able to read the BEDFORD files and extract the data so they 
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could be later on transferred to the Masscomp computer of the 
[RoeD. 
2. Code Analysis 
A program written in RM/ FORTRAN can read a data file 


written in one of the following forms: 


sa FORMATTED files consisting entirely of formatted 
records. 

2. UNFORMATTED files consisting entirely of unformatted 
records. 


The way these files are accessed/addressed determines how the 
records are read or written from a file. Furthermore, 
depending on the order in which the data are arranged inside 
the files, we can have sequential or direct access. The 
sequential files contain no logical gaps, can vary in length 
and can contain an End of File record. The direct access 
files can either be written or read in any order since they 
are random access files and they have a record number 
identifying the logical position of the record in the file. 

Our first attempt was to read the BEDFORD files using the 
sequential access method. The reason we did this was that by 
using a utility program (Norton Utilities) we recognized End 
of Files markers between the data (ASCII characters 1A and 
1D) which looked like Figure 9.1. 

Our attempt was unsuccessful due to the fact that, 
although the beginning headers are formatted in sequential 


records of 128 bytes each, the following data are records not 
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having a specific length but rather varying in length for 
every line representing a separate detector as described in 
Chapter V. 

The formatted method caused the integer data to be read 
as formatted characters instead of integers. For that reason 


every time an attempt to compile the program was made we 







encountered an error "-~1" meaning that an "End of File” was 
encountered on the reading attempt. 
Record n 


i 


Figure 9.1 Formatted Sequential Data 








We changed the method of accessing the data to the 
Unformatted Sequential method with unsuccessful results once 
more due to the fact that the data records this time were not 
sequential. 

We then tried the Unformatted Direct Access method and 
once more we had unsuccessful results. But this time there 
was no error in the method but we had an oversight error. 
When we opened the file as a direct formatted access record 
we should assign a record length defined by the “RECL" 
specifier in the open statement. The "“RECL" specifier was 
first set to one so we would be able to read the header 


characters one at a time. However the data matrix following 
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waS in integer*2 format and trying to read that format one 
byte at a time made the attempt unsuccessful. 

The problem was solved by opening the file once by using 
record length equal to one (RECL=1) in order to read the 
headers and then, after closing the file, re-opening it again 
using a record length equal to two (RECL=2). The reason we 
did this was to bypass the problem of the integer*2 of the 
data so the computer could read them. The file was opened by 
uSing as a dummy argument the file’s name (BEDFORDI, 
BEDFORD2, BEDFORD3). The format was “Unformatted”, the status 
of the file was "OLD" and the file access was “DIRECT”. More 
information about the code is presented in Appendix B. 

LSB MSB 
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Figure 9.2 DOS data SWAP for integer*2 data 
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Figure 9.3 BMAP data SWAP for integer*2 data 


We then had another problem to solve. The data were 


encrypted in the DOS basic coding for the integer*2 format 
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presented in Figure 9.2 with a small alteration. That 
alteration was that the most significant byte and the least 
significant byte were interchanged as presented in Figure 
oe. 

The Fortran program in Appendix B uses a subroutine 
called "SWAP" to achieve this. What that subroutine does is 
to reverse the order of bytes and replace the least 
Significant byte with the most significant byte. 

Finally after this problem was solved the data were read 
successfully and were transferred into three output files 
called “BEDFORD1.OUT", “BEDFORD2.OUT" and “BEDFORD3.OUT. 

We then used a normal modem communications program 
between computers called "“PCTALK 4" to transfer the three 


files to the Masscomp computer of the NPS-IRSTD system. 
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x. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

The BMAP program is a program with analysis, modeling and 
plotting capabilities of various infrared background scenes. 
It is created in support of the NSWC clutter computer code 
Originally designed for the Navy’s fleet defense systems and 
especially for analyzing a target’s background signature 
against a clutter background. 

The program was initially provided with three background 
infrared scenes called BEDFORD 1, 2, 3, which are typical 
representatives of the NATO Standard Format of recording and 
storing infrared scenes. 

We evaluated the BMAP program to find out its 
capabilities. The program proved to be excellent in analyzing 
and modeling the infrared data of the provided files with the 
program. For the analysis it used 2D Fast Fourier Transform 
comparison and other processing techniques, a version of the 
LOWTRAN propagation code and high resolution graphics for the 
screen including excellent plotting in a normal printer. We 
then evaluated the AGA DISCO program of the IF 800 computer 
and discovered a subroutine called "“REDA” which can be used 
to transfer the AGA infrared image file data to the IF 800 
memory and then by writing another computer code in BASIC or 


FORTRAN to transfer the same data to the IBM AT computer. 
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We evaluated the present data file form of the NPS/ IRSTD 
MASSCOMP computer for use of the BMAP program with a future 
application of transferring the data from the AGA through the 
IBM AT computer to the masscomp. 

We then attempted to transfer the three data sets from 
the BMAP program to tne Masscomp computer. It must be 
mentioned that the BMAP program is not limited to only the 
three BEDFORD data sets. It can also be used with all the 
data sets, now currently available from the Naval Research 
Laboratories, describing various infrared atmospheric scenes. 
Since these files are recorded and stored in the same method 
and format as the BEDFORD files it is assumed that the BMAP 
program will be able to perform its powerful capabilities 
with these files too. 

These data sets now can be transferred to the Masscomp 
computer of the IRSTD for further processing in conjunction 
with incoming data from the IRSTD sensor now on its final 
engineering development. In such away we will be able to 
combine the infrared target picture of the IRSTD with one of 
the infrared background scenes that will be stored inside the 
Masscomp computer of the IRSTD. 

For being able to do this we had first to transfer the 
BEDFORD files (typical Nato Standard Format infrared scenes) 
from the compatible IBM AT computer to the Masscomp computer. 
These files were first coded in a format described in the 


previous chapters which did not permit their transfer in an 
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understandable form. A computer algorithm had to be made that 
would be able to decode the data in to a simple alpha- 
numerical format that would be able to be transferred to the 
Masscomp computer with a modem. 

We then developed a computer software program in 
RM/FORTRAN, as described in Chapter IX and Appendix B, which 
was able to decode the BEDFORD files and rewrite them into 
three output files in alpha-numerical form. 

The output files created under this program were then 
transferred to the Masscomp computer using a modem and a 
Simple software transferring program called "PC-TALK 4” for 
internal storage and further processing with the geia ode 


infrared data. 


B. RECOMMENDATIONS 

Since we know that the Naval Research Laboratories have a 
great variety of infrared background scenes written according 
to the BMAP/BEDFORD format we recommend that these files be 
obtained for further analysis using the BMAP program. In 
that way we will not only explore the NRL computer library of 
infrared scenes but.will also create one of our own with the 
same data on the IBM compatible computer. 

It is then recommended that these infrared scenes be 
transferred to the Masscomp computer of the IRSTD, using the 
same developed RM/FORTRAN algorithm for further processing 


with the IRSTD data. 
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A further development of this thesis will be the transfer 
of infrared scenes from the small portable AGA Thermovision 
780 system to the Masscomp computer. To do this it is 
recommended that a computer code be written in BASIC or 
FORTRAN calling the "“REDA" subroutine in the IF 800 computer 
to read the infrared scenes of the AGA (already in storage in 
diskettes) in to memory. They can then be transfered to the 
IBM compatible AT computer using the "REDA" subroutine, as 
described in Chapter VII, and then, using the RM/FORTRAN 
algorithm, transfered to the Masscomp computer for further 
processing 

By doing this we should be able to make a comparison of 
the infrared scenes obtained with the AGA and with the IRSTD 
to make a rough estimation of what the IRSTD is really seeing 


and improving the capabilities of both systems 
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APPENDIX A 


LOWTRAN (LOWIN-LOWPLOT) INPUT FILES 


Model Atmosphere 

Type of Atmospheric Path 

Mode of Execution 

Temperature & Pressure Altitude Profile 
Water Vapor Altitude Profile 

Ozone Altitude Profile 

Radiosonde Data are to be Input 

Suppress Profile Output 

Temp at Boundary (.000 - T @ ist level) 


Surface Albedo (.000 - Blackbody) 


Extinction Type & Default Range 
Seasonal Arosol Profile 

Aerosol Profile & Extinction Type 
Air Mass Character 

Inclusion of Cirrus Attenuation 
Use of Army (VSA) for Aerosols 
Surface Range (.000 - Default) 
Current Wind Speed (m/sec ) 

24-Hr average Wind Speed (m/sec) 


Rain Rate (mm/hr) 


Meteorologic Data 
Horizontal Path 


Transmittance 


No 
No 
-000 


~000 


No Aersol Attenuation 
Determined by Model 


Background Stratosphric 


0 
No 
No 
~000 
~000 
-000 
~900 


Plot Type 


ny pe of X Axis 

Type of Y Axis 

Number of Decimal Digits for Y Axis 
Length of X Axis (in inches) 
Beginning Wavenumber/Waveiength 
Ending Wavenumber/Wavelength 

X Axis Annotation Interval 
Length of Y Axis (1in inches) 
Minimum Transmittance/Radiance 
{ 

‘Maximum Transmittance/Radiance 


Y Axis Annotation Interval 


Altitude (km) 

Pressure (mb) 

Ambient Temperature (degrees) 
Dewpoint Temperature (degrees) 
Relative Humidity (7%) 

Water Vapor Density (gm/m3) 
Ozone Density (gm/m3) 

Path Length (km) 

Initial Frequency (wavenumber ) 
Final Frequency (wavenumber ) 


Frequency Increment (wavenumber ) 


en Sp See 
Linear 
Linear 
i 
7.0000 
2000.0000 
4000.0000 
400.0000 
6.0000 
- OOE +00 


. OOE+00 


-0O0E+00 


-000 
-000 
-000 
- 000 
~000 
-000 
» 000 
2000 
2000.000 
4000 .000 


23-000 


APPENDIX B 


RM/FORTRAN COMPUTER CODE TO EXTRACT THE DATA FROM THE BEDFORD 
FILES OF THE BMAP SOURCE CODE 


Cc this program reads files partially successfully if they 


C are at least 30 lines long 


* 


this subroutine reads a matrix 


* 


RFILE has been changed to read binary rather than formatted 


K variable names 


*K FILENAME name of data file 

* X lst variable 

* Y 2nd variable 

* Z 3rd variable 

* NPTS number of data points 
*k inputs 


character filename*8, var*l,array(256)*l1 
* outputs 
real *4 x( 6000) 
* local 
logical here,there, everywhere 
integer Inunit, errorfile, jg kinpeseival (10. 
integer*2 data,cdata, bpdata(5936) 


character infile*8, dummy*10 


eS 


print *,’Input filename to open ’ 
read(6,’(a8)’) filename 
* 
* Assign the unit number 
famnit = J 
movetormat(i8) 
Cc print *,’input record length’ 
s=Q 
j=1 


3 


Pernt *,’ j= a 
open(unit=2,file=’ bedford. out’ ) 
* open file 
infile = filename 
open Grae = aus: 
file=filename, 
lostat=errorfile, 
status=’old’, 
form=’unformatted’, 
access=’direct’, 
recl=)) 
if (errorfile.ne.0) then 
print *, ’Error opening ’,filename 
close (Inunit) 
zo to 1000 
end if 


inquire(file-filename, 


as 


1OStav=1os, 
exist=here, 
opened=-there, 
access=dummy ) 
print *,’filename,ios,exist,opened, access’ 


print *,filename,ios,here,there,’ > , dummy 


x read heading 


38 


200 


on 


continue 

k=256 

print *,’input value of record length of header’ 
read (6,100) k 


3 


print *,’read ival from Inunit 

Print: sa 

doi { j= ier 

read(inunit,rec=j) var 

array(j)=var 

write (6,200) var 
write(2,200) var 

format(\al) 

continue 

print *,’ rewrite’ 

write (6,*) array 


close(Inunit) 
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c read data 


Cc print *,’input record length of data-’ 
Seep rint *, ° (i.e., number of bytes per data point):’ 
S read (6,100) 3 
j=2 
femmes. Ja  *., 3 
* open file 
infile = filename 
open Guam = Fnurnit, 


file=filename, 
iostatz=errorfile, 
status=’old’, 
form=’unformatted’ , 
access=’direct’ , 
recl=}3) 
if (errorfile.ne.0) then 
print *, ’Error opening ’,filename 
close (Inunit) 
go to 1000 
end if 
print *,’input starting data address’ 
read (6,100) j 
Prinwer, | ~ 
demoel=1, 5936 
read(inunit, rec=1+128) data 


call swap(data,cdata) 


ily 


300 


1000 


if (edata, 2 lt. hme -s 
bpdata(l)=cdata 


continue 


print data 

do 6 “L=O07 e900, iL 

write(6,300) bpdata(itl),bpdata(itZ)jbpearalia, 
bpdata(it4), bpdata(it+5), bpdata(i+ts6), bedatali+/e 
bpdata(it+8), bpdata(it+9), bpdata(itiO), dbpdata(itil) 
,bpdata(it+l2), bpdata(it13),bpdata(i+1l4),bpdata(i+15) 
write(Z,300) bpdata(itl),; dpdata(a-2)4 bedatal(iza @ 
bpdata(i+4),bpdata(it+5),bpdata(it+6),bpdata(it7), 
bpdata(it8),bpdata(it+9), bpdata(itidQ), bpdata(itil) 
continue 
Lormav( ix, bolo. ix) 


print *, s,’ negative numbers in data’ 


CleseeGlnunat) 

stop 

end 

subroutine swap(b,c) 
integer*2 b,c 

C=tsne (bh, 6) 
recur 


end 


APPENDIX C 
KEY FUNCTIONS FOR THE BMAP PROGRAM 
DISPLAY MODULE 


F2 key CROOK ICRI OK OK OK OK K KOK 


F2 key 


F3 key 


F4 key 


mo key 


F6 key 


Brings the DISPLAY MENU model 
and: 

Data Select 

Cursor 

Threshold Image 

Zee Let 

Ber Ot 

un 

Exits the Program 


OKO OKC OK OOOO COOK OOK KKK 
Gets a new set 
of images. 


FORO OOK OOOO OOO OR OIOK OK KOK 2K 
Brings the ANALYSIS MENU model 
and: 

Histogram 

Edge Operators 

Image Metrics 

2D PSD 

Filter PSD 

Drenal Frroc. 

Exits the Program 


OOOO CK COO OOK OOK IK OK KOK 
Points and 

reports pixel 

values. 

Terminates 

With the ESC key. 


OOK CORO OK RK OK KKK 
Brings the MODELLING MENU model 
and: 

Monte Carlo 

Structure Model 

LOWIN 

LOWTRAN 6.1 

LOWPLT 

DIR 

Exits the Program 


FORO OK KOCK OK IOI OK OKO OK CK OK KK KKK 
Changes the thresholds 


Zo 


F7 key 


F8 key 


FS key 


FO key 


for the imace caiaseray, 
OOO OOOO OOK KOK OKO OK KK CK kK 


2D Plotting 


FOOCIOOOOOOR OOOOC OOICOOOR I O ROK 
SUP ELotling 


ARC RCICRIC IOI KOIKE OK IC HOOK KK OC 


Allows the 

examination of 

a directory 

during program 

execution. 
OOOO OKO OK KIO KR OK OI OK kok ak 
DEOPS. brorccamn 

execution and 

returns to DOS. 
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Fi key 


F2 key 


F3 key 


F¥4 key 


F5 key 


F6 key 


MODEL MODULE 


SOCIO CIO OK OK KKK 
Brings the DISPLAY MENU module 
and: 

Data Select 

Cursor 

Threshold Image 

20) ieee 

SUE eloc 

iar 

Exits the Program 


OKC OCC HOR OO OK IO OK OK OK 
Monte Carlo 


ocattering 
Model 


OOK ORO OK OOO CKO OK OKO OK KOK 
Brings the ANALYSIS MENU module 
and: 

Histogram 

Edge Operators 

Image Metrics 

Zier oD 

Filter PSD 

Signal Proc. 

Exits the Program 


ORCC COO OOOOK IORI KK 
Cloud Structure Model 
--NOT CURRENTLY IMPLEMENTED 


COCK CCCI CIO ICR CORIO IOI OK KK 
Brings the MODELING MENU module 
and: 

Monte Carlo 

otructure Model 

LOWIN 

LOWTRAN 6.1 

LOWPLT 

DIR 

Exits the Program 

and returns to DOS 


OCC AO IC OK kK aC KOK 
LOWIN 


See the PC-TRAN manual [{Ref.23] 
and the Lowtran-6 manual [Ref.15] 


eA 


F7 key CCC OOO OOK OK IO IOK IC 
LOWTRAN 6.1 


See the PC-TRAN manual [Ref.23] 
and Lowtran-6 Manual [Ref.15] 


F8 key COO OK OK KR OK 
LOWPLT 


See the PC-TRAN manual [Ref.23] 
and Lowtran-6 Manual [{Ref.15] 


FS key CCK COKE KCK 2K 
Allows the examination of a 
directory during program execution 


FO key [OCCA OK IO KKK 
Stops program execution 
and returns to DOS. 
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Fi key 


F2 key 


Ho key 


F4 key 


WS key 


Bo, Key 


ANALYSIS MODULE 


FORO IOK OK OK OK 
Brings the DISPLAY MENU model 
and: 

Data Select 

Cursor 

Threshold Image 

Zp Plot 

5 Diao re 

ae 

Exit Program 


SOCIO CK ICRC ICO OOOO IK AOK 400K 
Histograms the 
cloud image. 


OKO OOOO OOK OOK KOK OK 
Brings the ANALYSIS MENU model 
and: 

Histogram 

Edge Operators 

Image Metrics 

2DaeoD 

Hanecer Pol 

Saenali Proc. 

Et r roeram 


FOCI OK OK aK 
Calculates 

metrics for the 

cloud image. 


NOT CURRENTLY 
IMPLEMENTED ! 


FOCI IOC OIC KOK OK 
Brings the MODELLING MENU 


Monte Carlo 
Structure Model 
LOWIN 

LOWTRAN 6.1 
LOWPLT 

DIR 
RXiterrogram 


SOCIO CCR OK OK HOCK ROK CO OK ORO HOI K 
Edge Operators 


NOT CORRENTLY 
IMPLEMENTED ! 


ae 


F7 key 


F8 key 


FS key 


FO key 


OOOO COICO OOK OR OOK KK 
Calculates and 

displays the 

2D Power 

Spectral density 


OOO OOOO OOK OKO OOK IK 
Performs 2D 

filteringwon 

the cloud image 

and displays 

the results. 


FOO CICK OOOO IOI OOK OIOKOK 
Signal Process 


NOT CURRENTLY 
IMPLEMENTED ! 


OOOO OOO OOO OOOO OK OKOK 
Stops program 

execution and 

returns tors. 


Wate 


APPENDIX D 


COMPUTER CODE SUPPLIED BY ONTAR COMPANY TO READ THE BEDFORD 
DATA FILES FROM THE BMAP CODE 


integer*2 data(371,16),i2 
characterx*l head(128),filenm(65),c1(2),chsv 
integer*2 openfh,rdmem,closfh,ier,nfh 


equivalence (c1,i12) 


a write(6,*) ’ Input Filename: ’ 
call gasciz(filenm) 
ier=openfh(0, filenm, nfh) 
if(ier .eq. 0) go to 2 
write(6, x) File error #’, ier 
Zomaro. 1 


Z write(6,*) ’ ’ 


Cc 
Cc read image 
Cc read two headers 


O 


ibuf=iaddrs (head) 

tfer=rdmem(nfh, ibuf,128) 

if(ter .ne. 0) go to 798 

write(*,30) (head(i) , i=1,128) 
30 format (1x,64a1/1x,64al) 


ier=rdmem(nfh, ibuf,128) 


AS 


101 


100 


O 


200 


if(ier .ne. 0) go to 798 


write(*,30) (head(i) , 
ibuf-=iaddrs(data) 

read data 

GomlOO yates 


ier=rdmem(nfh, ibuf,742) 


Law hZ8) 


iffier .ne. 0) wo to 798 


byte reverse 


Go VOd ses =1e53 7el 
i2Z=data(j,i) 
chsv=cl1(1) 
e1 (jel (2) 
Cl (2 ieensy, 
data(j,1i)=1i2 
continue 
ibuf=ibuf+742 


continue 


write data to screen 


write(6,*) ’ Input channel # 1-16, 


read (5,*) ichan 
iene Tena 


tf. wed. 0). isveo 
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O: “to> Hecate 


3 


738 


0 


1000 


100 


ioe enor leat, Lom goto 200 
write(6,*) ’ Channel # ae 
write(6,*) (data(j,i),j=1,371) 
WELte(6,*) * 

Poto 200) 


Stop 


write(6,*) *® read error #’ ,ier 


stop 


end 


subroutine gasciz(name) 


character*l name(65) 


Head os, lUOU =i mame(i) , i=1,65) 
format(65al1) 

do 100 i=1,64 

ifcnametin .eq. «= ) go to 1 
eConcenue 

i1=65 

name(i)=char(0) 

meter 


end 
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